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ABSTRACT 

This  mathematical  model  is  designed  to  provide  predicted  dynamic  per¬ 
formance  data  of  shorebased  rotary-hydraulic-type  aircraft  arresting  gears. 
A  Navy  Kodel  E-28  arresting  gear  is  used  for  specific  comparison  between 
predicted  results  of  the  computer  solution  and  actual  test  results.  The 
siirulation  of  an  arrestment  of  a  vehicle  under  a  particular  set  of  condi¬ 
tions  is  accomplished  by  putting  information  (data)  into  the  computer. 

The  input  data  specifies  values  for  the  installation  geometry  and  mechan¬ 
ical  properties  of  the  arresting  system  and  the  test  vehicle.  Predicted 
■iyrsmi-  "alues  of  forces  and  motions  of  the  test  vehicle,  purchase  system, 
arc’  t*pe  reel  are  printed  out  versus  time  at  a  predetermined  incremental 
time. 

Inis  report  is  a  phase  report  on  tie  development  of  the  model  and 
contsir?  the  early  analytical  design  approaches,  the  most  current  analyt¬ 
ic  -1  reproach  with  the  computer  program,  and  instructions  for  execution 
of  the  computer  program. 
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I  INTRODUCTION 

A.  The  E-28  Mathematical  Model  was  designed  and  is  being  developed 
by  the  Computer  Division,  Engineering  Department,  NATF,  to  simulate  the 
performance  characteristics  of  the  E-28  arresting  gear  and  any  similar 
shorebased  rotary -hydraulic  arresting  gear,  such  as  the  BAK-13  or  44B-2D. 

B.  The  ultimate  objectives  of  the  design  are: 

1.  To  perform  parameter  predictions  which  aid  in  arresting-gear 
component  design,  equipment  modification,  and  the  determination  of  gear 
performance  changes  due  to  various  proposed  and/or  actual  installation 
configurations. 

2.  To  reduce  the  test  time  of  a  calibration  test  program  by 
generating  performance  data  which  could  act  as  "fill  in"  data  once  the 
upper  and  lower  limits  of  the  engagingrspeed  and  vehicle  gross-weight 
spectrums  have  been  established  through  actual  arrestments,  provided 
good  agreement  exists. 

C.  This  report  has  been  prepared  to  document  the  model  in  its  cur¬ 
rent  state  of  development  and  to  discuss  the  additional  capability  neces¬ 
sary  to  achieve  the  above  objectives  (paragraphs  B1  and  B2) . 
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II  DESCRIPTION  OF  THE  E-28  ROTARY -HYDRAULIC -TYPE  ARRESTING  GEAR 

A.  Purpose:  The  Navy  E-28  arresting  gear  is  a  shorebased  emergency 
arresting  gear  designed  to  arrest  all  U,S.  Navy  arresting -hook-equipped 
airplanes  under  conditions  of  aborted  takeoff  or  landing  overrun. 

•  Capabilities:  An  airplane  engaging  the  arresting  gear  will  be 
stopped  within  a  runout  distance  of  approximately  1,000  feet.  The  maxi¬ 
mum  energy  absorbing  capacity  of  the  gear  is  76  million  foot-pounds 
(nominal).  Engagements  can  be  made  from  either  runway  direction  and  at 
points  up  to  40  feet  on  either  side  of  the  runway  centerline. 

*  * 

C.  General  Description  of  the  Arresting-Gear  Operation:  As  shown 
in  Figure  1,  two  identical  energy  absorber  units  and  runway-edge  sheaves 
are  ltKated  on  opposite  sides  of  the  runway  and  are  connected,  through 
nylon  purchase  tapes,  to  a  steel  wire-rope  deck  pendant.  Arrestment  of 
a  landing  aircraft  is  accomplished  by  engagement  of  the  aircraft  arrest¬ 
ing  hook  with  a  pendant  stretched  across  the  runway.  The  attached  pur¬ 
chase  tapes  are  pulled  off  the  six-foot-diameter  drum  on  each  arresting 
gear.  Each  drum  is  splined  to  a  shaft  which  turns  a  vaned  rotor  between 
vaned  stators  in  a  housing  filled  with  a  water/glycol  mixture.  The  tur¬ 
bulent  fluid  resistance  caused  by  the  stator  and  rotor  interaction  (water 
brake)  decreases  the  rotational  speed  of  the  drums,  thereby  slowing  down 
the  purchase-tape  payout  which  in  turn  applies  a  braking  force  on  the 
aircraft.  The  ensuing  fluid  turbulence  converts  the  landing  aircraft's 
kinetic  energy  into  heat.  A  codling  system  is  provided  to  dissipate  this 
heat  during  rapid-cycle  operations.  After  the  aircraft  has  been  safely 
brought  to  a  stop,  and  the  arresting  hook  disengaged,  the  pendarit  and  nylon 
tapes  are  returned  to  battery  position  by  an  air-cooled  gasoline  engine 
driven  retraction  system. 
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III  BACKGROUND  OF  MODEL  DEVELOPMENT 

A.  General:  The  development  of  the  model  from  the  initial  working 
model  to  the  current  model  has  mainly  involved  modifications  of  the  pen¬ 
dant/tape  geometric  configuration  that  occurs  during  an  arrestment  and  is 
used  as  the  basis  for  parameter  generation. 

B.  Design  Progress  to  Date 

1.  Initial  Design:  The  initial  E-28  mathematical  model  produced 
dynamics  of  the  gear  directly  from  the  motion  of  the  arresting  aircraft. 

The  geometry  of  the  pendant/tape  pattern  was  assumed  to  be  triangular  in 
shape.  The  pendant/tape  was  assumed  to  be  in  a  straight  line  proceeding 
from  the  arresting  hook  to  the  runway -edge  sheaves  throughout  the  entire 
arrestment.  The  early  results  of  model  runs  indicated  that  a  more  sophis¬ 
ticated  approach  was  necessary. 

2.  Interim  Design 

a.  The  kink  wave  is  a  triangle-shaped  deformation  of  the  run¬ 
way  pendant  that  is  generated  from  the  impact  of  the  arresting  hook  and 
represents  the  motion  of  stress  propagation  in  the  pendant/tape  setup. 

b.  The  idea  of  designing  kink -wave  phenomenon  into  the  math 
model  originated  after  studying  reports  written  by  F.O.  Ringleb  (refer¬ 
ences  (a)  and  (b))  concerning  cable  dynamics.  The  addition  of  kink -wave 
motion  td  the  program  design  was  necessary  to  more  accurately  predict 
E-28  arresting-gear  performance,  especially  arresting-hook-load  and  tape- 
tension  values,  for  the  initial  part  of  an  arrestment.  The  introduction 
of  kink-wave  geometry  into  the  model  forced  the  dynamics  of  the  arresting 
system  to  be  generated  with  respect  to  the  motion  of  the  kink -wave.  The 
various  methods  chat  were  devised  to  simulate  actual  kink-wave  paths 
which  hold  for  the  entire  simulation  of  an  arrestment  are: 

(1)  arresting -hook-point-to-sheave  motion 

(2)  arresting -hook-point-to-sheave/ sheave-to-arresting  - 
hook-point  motion 

(3)  arresting-hook-point-to-sheave/sheave-to-arresting- 
hook-point/arresting-hook-point-to -sheave  motion 

3.  Most  Recent  Design:  The  best  results  that  have  been  obtained 
to  data  are  from  the  current  model.  This  approach  assumes  that  the  kink 
A-ave  travels  from  arresting -hook  point  to  sheave  repeatedly.  That  is,  when 
ct-  t  program  determines  that  the  kink  wave  has  reached  the  runway-edge  sheave 
ar.d  l  *  reaay  to  "bounce  back"  toward  the  hook,  program  logic  forces  a  new 
kiOk  .,s;s  to  emanate  from  the  hook  and  progress  coward  the  sheave.  Although 
not  t n  agreement  with  the  classical  kink -wave  motion  described  by  Ringleb  m 
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reference  (a),  this  method  has  proven  to  give  the  most  accurate  predic 
tion  results  to  date.  The  results  of  this  method  are  compared  to  the 
results  of  earlier  design  approaches  in  the  form  of  a  composite  plot  of 
arresting-hook  load  versus  time  in  Figure  2. 
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IV  ACTUAL  ARRESTMENT  DESCRIPTION 

A.  The  purpose  of  the  arresting  system  is  to  dissipate  the  engaging 
kinetic  energy  of  an  aircraft.  This  is  accomplished  by  transmitting  the 
vehicle's  energy  through  the  purchase  tapes  into  the  tape-reel  absorber 
units  where  the  energy  is  converted  into  heat  by  the  action  of  fluid  tur¬ 
bulence.  During  the  course  of  an  arrestment  there  are,  therefore,  three 
main  types  of  motion  caused  by  the  action  and  reaction  of  the  system 
components : 

1.  Tape-reel  and  rotor  motion 

2.  Pendant/ tape  (kink-wave)  motion 

3.  Vehicle  motion 

The  force  interaction  between  these  three  types  of  motion  is  the  principal 
basis  for  the  mathematical  model  simulation  of  an  arrestment. 


B.  The  following  is  an  account  of  how  each  type  of  motion  changes 
as  the  arrestment  proceeds: 


1. 

Before  Impact 

a . 

Arresting  Gear  -  motionless 

b. 

Pendant/Tape  -  motionless 

and  weight 

c. 

Vehicle  -  approaching  pendant  at  a  predetermined  speed 

2. 

At  Impact 

. 

Arresting  Gear  -  motionless 

b. 

Pendant/Tape  -  set  in  motion  by  action  of  kink-wave 

generated  upon  impact 

c.  Vehicle  -  speed  at  impact  is  the  engaging  speed  for  the 

arrestment 

3.  After  Impact  (period  of  arrestment):  The  initial  braking 
force  on  the  vehicle  is  slight  and  is  imposed  by  tensions  due  to  purchase- 
tape  elongation.  Once  the  tape  reels  are  in  motion,  however,  the  main 
retarding  (braking)  force  is  due  to  the  hydraulic  brake  connected  to  the 
tape  r*>ol .  A  typical  arrestment,  therefore,  will  be  described  in  two 
sections:  the  first  called  the  dynamic  region  where  there  is  only  slight 

retard, if  inn  of  tin*  vehicle  due  to  kink-wave  motion  causing  purchase-tape 
elongation;  and  the  second  called  the  hydraulic  region,  where  the  main 
retarding  on  tho  vehicle  are  encountered. 
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a.  Dynamic  Region 

(1)  Arresting  Gear  •  tape  reels  start  to  rotate  due  to 
tension  instilled  by  kink  waves.  Fluid  pressure  in  absorber  is  still 
not  affected. 


(2)  Pendant/Tape  -  kink  wave  is  reflected  off  of  runway* 
edge  sheaves  and  travels  back  and  forth  repeatedly  from  arresting-hook 
point  to  3heave,  causing  the  kink-wave  "humps"  in  the  recorded  tape- 
tension  values. 

(3)  Vehicle  -  speed  starting  to  feel  effects  of  tension 

puli  of  tapes. 

b.  Hydraulic  Region 

(1)  Arresting  Gear  -  tape  reels  are  increasing  speed 
rapidly  because  they  are  feeling  directly  the  puli  of  the  vehicle's 
weight.  Fluid  pressure  in  Absorber  unit  Increases  steadily, 

(2)  Pendant /'Tape  •  kink  waves  have  more  or  less  damped 
out  although  effects  on  tape  tension  are  still  noticeable.  Tape  tensions 
are  due  almost  entirely  to  the  retarding  force  the  taps  reel  fluid  has 

on  the  vehicle  via  the  purchase  tape. 

4.  End  of  Arrestment:  When  ail  of  the  vehicle's  kinetic  energy 
has  been  transformed  into  heat  in  the  absorber  units  through  the  action 
of  fluid  turbulence,  the  vehicle  comes  to  a  stop  and  the  arrestment  is 
complete.  Reel  acceleration,  arresting-hook  load,  absorber  pressure, 
and  tape  tensions  reached  their  maximum  values  in  the  hydraulic  region 
of  the  arrestment  and  then  decreased  to  their  initial  pre-impact  values*-. 


^Thi?  is  generally  true  although  sometimes  the  maximum  values  of 
arresting -hook  load  and  tape  tensions  can  reach  their  peaks  in  the  dynamic 
region  for  certain  arrestment  runway  configurations.  Also,  due  to  tape 
stretch  there  may  be  a  lead  in  the  system  when  the  vehicle  is  stopped 
announcing  for  small  arresting -hook-load  and  tape-tension  values. 
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V  SIMULATED  ARRESTMENT  DESCRIPTION 

A,  General;  The  computer  program  developed  is  a  finite  incremental 
analysis  of  arresting-gear  dynamics.  The  simulation  of  a  test  event 
(aircraft  arrestment)  is  accomplished  by  relating  mathematically  the 
interreactions  between  vehicle,  kink-wave,  and  tape-reel  motions.  Param¬ 
eter  computation  sections  of  kink-wave  velocity,  tape  tension  and  elon¬ 
gation,  tape-reel  acceleration,  etc.,  which  have  been  developed  in 
separate  studies  are  applied  collectively,  while  conforming  to  the  pre¬ 
determined  shape  or  pattern  of  kink-wave  design,  to  obtain  performance 
results. 

B.  Calculation  Section  Logic 

1.  The  calculation  section  logic  consists  of  equation  segments 
which  compute  specific  physical  properties  of  the  arrestment  such  as 
forces,  torques,  motions,  and  distances.  Each  segment  was  developed 
independently  and  then  combined  to  provide  for  the  calculation  of  the 
key  motion  generation  parameters--vehicle  and  reel  accelerations. 

2.  The  logic  that  begins  the  calculation  loop  and  generates 
simulated  motion  contains  equations  which  are  actually  simple  integration 
techniques.  Values  of  accelerations  (tape  reel  and  vehicle)  are  provided 
from  the  previous  pass  through  the  calculation  loop.  The  integration 
technique  is  applied  which  results  in  values  of  velocity  and  speed  (tape 
reel  and  vehicle),  and  Che  procedure  is  repeated  to  obtain  values  for 
tape-reel  revolution  and  vehicle  runout.  The  basic  loop  logic  is 
diagrammed  below: 


-VEHICLE  ACCEL 


START! 

OF 

ICALC 

R.OOP 


^VEHICLE 
ACCEL  ’ 


JAPE-REEL 

ACCEL 


^VEHICLE 

RUNOU'f' 


^  REEL 
ROTATION 


f 


TAPE  ON _ JTAPr _ ^  TAPE  ARRESTING 

DECK  ELONG  TI  E  SI  ON  *HQOK  LOAD 


-REEL  ACCEL 
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3.  An  explanation  illustrating  how  vehicle  speed  is  obtainid 
from  vehicle  acceleration  follows: 

letting  V2  *  vehicle  speed  at  end  of  time  increment  (VELNS2), 

Vt  ■■  vehicle  speed  at  start  of  time  increment  (VELNST), 

A2  *  vehicle  acceleration  at  end  of  previous  time  increment. 
(VACC2), 

A-  -  vehicle  acceleration  at  start  of  previous  time  increment 
(VACC), 

T  -  time, 

and  At  *  At. 

T 

Instantaneous  speed  can  be  approximated  by  a  value  of  average  speed  or 
V2  -  \  +  At  <At  +  A2). 

This  expression  represents  an  integrated  process  which  utilizes  the 
average  of  two  accelerations,  A^  and  A2,  within  a  specified  time  incre¬ 
ment,  AT,  to  obtain  an  average  speed  difference,  A V  ■  At  (A^  +  A2). 

4.  Ideally,  the  most  Accurate  simulation  can  be  obtained  when 
AT  approaches  zero,  then  AV  ■  df.  This  program  assumes  that  sufficient 
accuracy  is  obtained  by  letting  the  time  increment  At  ■  .001  second. 

5.  Passing  through  the  calculation  loop  is  then  a  matter' of 
jumping  from  one  logic  segment  to  the  next.  The  calculation  loop  con¬ 
tains  seven  main  logic  segments  which  compute  values  for: 

a.  Reel  and  vehicle  motion 

b.  Pendant  or  tape  kink -wave  velocity 

c.  Kink -wave  coordinates 

d.  Pendant/tape  tensions  and  elongations 

e.  Tape-reel  inertia  and  acceleration 

f.  Arresting -hook  load 

g.  Vehicle  acceleration 
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A  detailed  description  of  how  each  of  these  segments  is  designed  can  be 
found  in  Section  VI. 
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VI  PROGRAM  DESCRIPTION 

A.  General:  The  program  description  format  parallels  that  of  the 
program  listing  (Appendix  A),  The  description  is  divided  according  to 
the  principal  sections  of  program  legic.  Variable  titles  used  here  are 
identical  to  those  in  the  program.  A  complete  list  of  variable  names 
and  definitions  is  presented  in  Appendix  B. 

B.  Computer  Information 


1.  The  program  is  self-contained  and  therefore,  does  not  refer 
to  external  sources  for  information,  The  program  language  is  FORTRAN  IV, 
and  the  program  has  been  run  primarily  on  the  CDC  6600  computer,  however, 
it  can  also  be  run  on  an  IBM  360/65  computer.  The  following  is  run 
information  for  an  average  run  on  the  CDC  6600? 


Compilation  time 


6  seconds 


Run  time 

Memory  locations  used 


16  to  20  seconds  (with  .001  second 
calculation  time) 

10,000  words  (60  bits  per>word) 


Output  line  limit 


2,500  lines 


2,  Data  is  written  on  temporary  disc  storage  to  minimize  running 
time.  Normally,  a  run  is  made  on  priority  6  status  which  costs  33  cents 
per  system  second.  If  turnaround  time  is  not  a  critical  factor,  the  cost 
may  be  reduced  by  running  on  priority  0  at  20  cents  per  system  second. 


C.  Program  Sections 
1,  Input  Section 

a.  This  section  contains  a  dimension  statement  which  defines 
the  linear  array,  DATUM(32),  which  is  used  for  program  output  generation 
and  storage.  Every  parameter  value  of  output  data  is  allocated  a  space 
in  this  array  during  calculations.  Data  values  can  be  recorded  on  tape 
or  disc  storage  each  time  the  calculation  loop  is  completed  or  values  can 
be  recorded  at  a  chosen  time  interval.  Presently  the  recording  time 
interval  used  is  .01  second. 


b.  The  variable  JGEOR  is  initialized  to  zero.  JGEOR  is  a 
sequential  run  indicator  and  controls  the  number  of  program  runs  performed. 
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c.  Read  and  format  statements  of  input  data  complete  this 
section  and  are  separated  to  distinguish  arresting-gear  configuration 
data  from  event  information  data.  The  various  input  values  that  are 
read  into  the  program  are  listed  in  Figure  3,  which  is  a  sample  print¬ 
out  of  ail  input  values. 

2,  Initialization  Section 


a.  Initial  values,  that  is  values  of  variables  that  are 
required  for  the  first  pass  of  the  main  calculation  loop,  generally  fall 
into  three  main  catagories: 

(1)  Deck  jjeometry 

(2)  Mechancial  properties  of  arresting  gear  (reel,  tape, 

pendant) 

(3)  Vehicle  and  arresting-hook  properties 

b.  The  following  initialization  operations  are  performed  on 
variables  in  these  categories: 

(1)  Assignment  to  the  Calculation  Variable  Names  the 


Value  of  Zero 


(a)  All  variable  -names  that  before  pendant  impact 
have  the  value  of  zero--such  as  tape  elongation,  kink-wave  coordinates 
(angles  and  distances),  arresting-gear  motipn  variables,  deck-geometry 
values,  arresting-hook- load,  tape  tensions,  etc. --are  set  equal  to  zero. 

(b)  Other  variables  set  equal  to  zero  are  TIME,  the 
symmetry  flag  IFLGSM,  the  output  printer  line  counter  ICOUNT,  and  the 
linear  array  DATUM(32). 

(2)  Conversion  of  Units  of  Input  Variables:  To  simplify 
the  preparation  of  input  values,  their  units  are  those  of  standard,  com¬ 
monly  used  values.  However,  to  ease  calculations  through  conformity  of 
units,  some  input  values  must  be  converted  to  units  which  agree  with 
other  variable  units.  For  example,  vehicle  engaging  speed  VELENG  is 
normally  referred  to  in  knots  and  is  therefore,  input  this  way.  For 
calculations,  however,  speed  in  FT/SEC  would  be- more  suitable  and  there¬ 
fore  the  conversion  is  performed. 
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(3)  Calculation  oflntermediate  Input:  Parameter  vari¬ 
ables  that  ace  not  part  of  ihput  but  have  definite  values  aside  from  zero 
and  must  be  introduced  into  the  main  calculation  loop,  are  calculated 
from  input  data  or  companion  values.  An  example  is  the  initial  vehicle 
acceleration  VACC,  which  is  computed  from  the  vehicle  mass,  VMASS,  and 
forces  of  thrust  VTHKUS,  arresting -hook- load  VHOOK,  and  drag  VTDRAG. 

3.  Calculation  Section 

a.  The  purpose  of  this  section  is  to  compute  prime  parameter 
values  as  functions  of  time  for  output.  The  parameters  that  are  currently 
being  computed  are  listed  in  the  output  section,  page  16..  To  compute 
prime  parameter  values,  many  interrelated  equations  describing 

(1)  Vehicle  distances  and  motion, 

(2)  Tape  and  pendant  properties  and  geometry,  and. 

(3)  Tape-reel  rotation  and  motion 

are  drawn  upon  to  develop  changes  of  the  "key"  motion  generation  param¬ 
eters,  vehicle  and  reel  accelerations,  with  time.  These  motion  genera¬ 
tion  parameters  are  the  basis  of  the  calculation  loop. 

b.  Time  Generation:  The  impact  of  the  arresting -hook  with 
the  pendant  is  simulated  when  the  calculation  loop  is  entered  for  the 
first  time.  For  each  pass  through  the  calculation  loop,  TIME  (arrest¬ 
ment  time)  is  increased  by  one  time  calculation  Increment. 


c.  Port  and  Starboard  Calculations:  Port  and  star-board 
parameter  values  can  be  calculated  separately  if  the  arrestment  is  OFF- 
CENTER.  If  the  symmetry  flag  indicates  an  ON-CENTER  arrestment,  calcula¬ 
tions  are  performed  for  port  parameters  only  and  set  equal  to  starboard 
values  for  each  pass  through  the  loop.  (See  Appendix  C,  Program  Flow 
Chart,  for  more  details.) 


d.  Main  Loop  Calculations 

(1)  The  main  calculation  loop  is  separated  into  five 

sections: 

(a)  Program  stop  logic, 

(b)  Port  parameter  calculations, 

(c)  Setting  of  port  parameter  values  to  starboard 
parameter  values  for  a  symmetrical  arrestment, 

(d)  Starboard  parameter  calculations,  and 

(e)  Final  bAlculations. 
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The  port  and  starboard  calculations  contain  the  3ame  equations,  differing 
only  by  the  variable  name  endings,  P  for  port  and  S  for  starboard.  There* 
fore,  only  the  port  parameter  calculations  will  be  described.  Variable 
names  that  depict  parameter  values  at  the  end  of  a  time  increment  end  in 
the  number  2. 


(2)  Program  Stop  Logic:  Two  separate  tests  are  performrd 
upon  entrance  into  the  main  calculation  loop  to  determine  if  the  program 
has  simulated  an  arrestment  completely  and  to  indicate  a  halt  to  program 
calculations: 


(a)  Port  or  starboard  two -block:  The  instantaneous 
radius  of  outer  wrap  of  the  purchase  tape  on  the  reel  ROWP  is  compared  to 
the  reel-hub  radius  RHUB  to  establish  if  the  tape  supply  on  the  reel  has 
been  depleted. 


(b)  Vehicle  Speed:  The  instantaneous  vehicle  speed 
VELNST  ,is  checked  to  see  if  it  has  reached  a  value  of  zero  which  would 
indicate  the  end  of  an  arrestment. 

(3)  Port  Parameter  Calculations 

(a)  Motion  calculations:  In  this  section,  motion  of 
the  vehicle  and  tape  reels  is  simulated,  Speed  and  velocity  (vehicle  - 
VELNS2,  reel  -  RSPP2)  are  obtained  by  averting  acceleration  values  of 
the  beginning  (vehicle  -  VACC,  reel  -  RACP)  a"hd  end  (vehicle  -  VACC2, 
reel  -  RACP2)  of  the  last  time  increment,  multiplying  by  the  time  incre¬ 
ments,  and  then  by  adding  this  value  to  the  speed/velocity  values  at  the 
end  of  the  previous  increment. (vehicle  -  VELNST,  reel.-  RSPP) .  Following 
this  same  integration  logic,  vehicle  runout  VRUN02;:.and  reel  revolutions 
RP0P2  are  obtained  from  speed/velocity  values.  The  radius  of  outer  wrap 
of  reel  tape  ROWP  is  calculated  from  the  amount  of  tape  that  has  left 
the  reel.  The  amount  of  tape  on  a  reel  at  any  time  RLTP  is  determined 

by  the  differences  of  the  total  length  of  tape  in  system  TOLENP  and  the 
amount  on  deck  TLENP  from  the  last  calculation  increment.  A  new  value 
for  tape  on  deck  TLENP2  to  account  for  the  motion  of  vehicle  and  reils 
is  then  computed.  See  Figure  4  for  description  of  variable  names  of  the 
runway  configuration. 

(b)  Kink-wave  location:  The  location  of  the  kink- 
wave  in  the  tape-pendant  configuration  is  established.  The  velocity  of 
the  kink  wave  is  a  function  of  the  media  in  which  it  is  traveling.  It 
is  necessary,  therefore,  to  determine  for  each  calculation  if  the  kink 
wave  is  traveling  through  the  nylon  tape  or  3teel  pendant  so  that  the 
proper  dynamic  equations  describing  its  motion  can  be  employed.  This 
test  is  accomplished  by  checking  the  three-dimensional  distance  HYPOTP 
which  is  the  distance  from  hook  to  kink  against  the  length  of  the  pendant 
on  the  port  side  of  hook  engagement  through  the  use  of  an  IF  statement. 
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(c)  Kink-wave  velocity:  The  kink-wave  velocity  for 
either  the  pendant  or  tape  is  calculated.  This  velocity  1*  computed  in 
terms  of  modulus  of  elasticity  (pendant  -  PENMOD,  tape  -  ETAPEP),  longi¬ 
tudinal  wave  velocity  (pendant  -  CPBND,  tape  -  CTAPEP),  and  the  trans¬ 
verse  impact  formula  approximation  (pendant  -  PEfSIG,  tape  -  TASIGP). 

(d)  Kink -wave  three-dimensional  coordinate  location: 

The  instantaneous  location  of  the  kink  wave  is  described  by  a  three-dimen¬ 
sional  coordinate  system  which  pinpoints  its  position  by  means  of  angles 
and  linear  distances  from  a  fixed  reference.  See  Figure  5  for  an  illus¬ 
trated  description;  of  this  coordinate  system.  Also,  the  instantaneous 
arresting -hook  elevation  HKELEV,  the  distance  of  the  kink  to  the  tail 
hook  HYPQTP,  and  the  distance  of  the* 'kink  to  the  runway-edge  sheave 
RESULP,  are  computed.  See  Figure  6  for  an  illustrated  description. 

(e)  Pendant  and  tape  elongation,  tension:  Total 
elongation  of  the  tape/pendant  configuration,  DELP,  is  computed  by  sub¬ 
tracting  the  length  of  the  tape/pendant  on  the  runway  of  the  previous 
time  increment,  TLENP2  and  PLENP,  from  the  new  length  of  the  tape/pendant 
on  the  runway  just  established  from  the  three-dimensional  distances, 

HYPOTP  and  RESITLP,  and  SPLITP.  See  Figures  4  and  6.  The  tape  elongation 
factor,  TELFP,  which  is  the  percent  tape  elongation  and  is  required  for 
substitution  into  the  tape  modulus  of  elasticity  equation  and  the  tape 
tension  equation  (see  pages  E-3  and  E-4,  Appendix  E  under  procedure),  is 
computed  by  subtracting  the  pendant  elongation  PHLP  from  the  total  elonga¬ 
tion  DELP,  and  dividing  by  the  product  of  the  length  of  tape  on  the  run¬ 
way  TLENP2  and  100.  Tape  tension,  DTENP,  is  computed  and  its  values  used 
to  calculate  a  new  value  of  pendant  elongation  PELP2,  which  is  a  function 
of  pendant  length  PLENP,  pendant  modulus  of  elasticity  PENMOD,  pendant 
cro3s-sectional  area  PENX,  and  tape  tension. 

(f)  Tape-reel  acceleration,  arresting -hook  load: 
Tape-reel  acceleration  RACP2,  is  determined  by  the  relationship  between 
the  torques  acting  on  the  reel  and  rotor  and  the  polar  moment  of  Inertia 
of  the  reel  and  rotor.  (See  Appendix  F  for  detailed  description  of 
equation  derivation.)  Arresting-hook  load,  H00KP,  is  established  by 
resolving  tape  tension  into  vector  components  in  the  direction  of  hook 
engagement,  (See  Figure  7,  page  31  for  diagram.) 

(g)  Test  to  determine  if  kink  wave  has  reached  the 
sheave  At  the  end  of  the  port  parameter  calculation  section,  a  test  is 
performed  to  determine  whether  the  kink  wave  has  reached  the  runway-edge 
sheave .  This  test  is  accomplished  by  comparing  XWAVEP,  which  is  the  dis¬ 
tance  from  the  centerline  of  engagement  to  the  kink  wave,  to  DHAP,  which 

i*  the  distance  from  the  centerline  of  engagement  to  the  runway-edge  sheave. 
When  the  test  indicates  that  these  values  are  equal,  the  kink  wave  has 
reached  the  sheave  and  the  three-dimensional  coordinates  and  angles  at 
j-j- 1 4  time  must  be  recorded  in  order  to  establish  a  new  path  for  the  next 
ki'V  wave  to  follow.  If  the  kink  wave  has  reached  the  sheave,  values  of 
CAMAXP ,  C4MA2P,  GAmAZP,  XWAVEP,  YWAVEP,  ZWAVEP,  and  HYPOTP  are  computed. 
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(4)  Arresting  Symmetry  Check:  If  the  symmetry  flag  IFLGSM 
indicates  an  ON-CENTER  arrestment  and  a  completely  symmetrical  .system 
installation  (IFLGSM  *  0),  port  parameter  values  are  set  equal  to  star¬ 
board  parameter  vaiues.  This  procedure  minimizes  computational  time  under 
computer  control  when  the  arrestment  is  symmefltlcal.  If  the  arrestment 
is  not  symmetrical  (IFLGSM  -  1),  the  computer  is  instructed  to  enter  the 
starboard  calculation  section.  This  test  is  performed  by  the  use  of  an 
IF  statement. 


(5)  Final  Calculations 

(a)  In  the  final  calculation  section,  VACC2,  the 
vehicle  acceleration  is  calculated  and  parameter  variable  values  that 
have  been  calculated  in  the  port  and  starboard  sections  (of  present  time 
increment)  are  stored  for  reuse  upon  the  re-entrance  of  the  main  calcula¬ 
tion  loop. 

(b)  Vehicle  acceleration  is  obtained  simply  by  apply¬ 
ing  Newton's  second  law  of  motion  and  using  the  fundamental  quation, 

Force  ■  mass  x  acceleration.  Vehicle  acceleration  is  equal  to  the  result¬ 
ant  forces  acting  on  it,  thrust  (VTHRUS),  arresting -hook  load  (VH00K), 
and  the  total  drag  force  (VTDRAG)  divided  by  its  mass  (VMASSO. 

(c)  Storage  of  parameter  values  for  use  in  the  next 
time  increment  is  accomplished  by  assigning  present  increment  variable 
values  which  are  suffixed  with  the  number  "2",  to  the  same  respective 
variable  names  not  followed  by  the  number  "2". 

(6)  Saving  of  Maximum  Values:  The  maximum  value  of  cer¬ 
tain  parameters  is  required  for  most  arresting  test  programs.  This  sec¬ 
tion  provides  a  means  of  obtaining  maximum  values  of  desired  parameters 
for  an  arrestment  along  with  the  times  that  the  maximum  values  occurred. 
Maximum  values  are  designated  by  the  variable  name  prefixed  by  the  letter 
"A".  The  following  maximum  values  are  currently  being  computed  for 
output: 


arresting-hook  load 

-  A7H00K 

tape  tension,  port 

-  ABTENP 

tape  tension,  stbd 

-  ADTENS 

vehicle  deceleration 

-  AVACCG 

reel  velocity,  port 

-  ARSPP 

reel  .velocity,  s.’  bd 

-  ARSPS 

reel  acceleration,  port 

-  ARAC? 

reel  acceleration,  stbd 

-  ARACS 

(7)  Main  Calculation  Loop  Output  Storage:  All  parameter 
values  that  are  to  be  output  are  stored  in  the  linear  array  DATUM(32) 
along  with  the  arrestment  time  TJME(DATUM  (9)),  and  are  recorded  on  tape. 
The  program  listing  statement  300  (see  Appendix  A),  shows  the  parameter 
names  along  with  their  assigned  array  position. 
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(8)  Run  Termination:  Run  termination  is  based  on  elapsed 
time  and  occurs  when  the  elapsed  time  (TIME)  exceeds  a  preset  arrestment 
time  limit. 

(9)  Kink -Wave  Coordinate  Generation:  The  last  section  of  < 
the  main  calculation  loop  generates  the  three-dimensional  kink-wave  coor¬ 
dinate  values  and  the  kink-wave  velocity  for  both  the  port  and  starboard 
kink  waves  along  with  the  corresponding  time  of  arrestment. 

4.  Output  Section 

a.  Program  output  consists  of  five  parts: 

(1)  Record  of  input, 

(2)  Kink -wave  coordinates, 

(3)  Maximum  values, 

(4)  Vehicle  motion  parameters,  arresting-hook  load  and 
tape  tensions,  and 

(5)  Port  tape  and  tape-reel  values, 

b.  The  record  of  input  data  is  output  on  printed  form  imme¬ 
diately  following  the  initialization  section  of  the  program.  See  Appen¬ 
dix  A  for  a  sample  output  listing. 

c.  Kink -wave  coordinates  are  current  values  that  are  printed 
after  every  calculation  iteration  until  the  run  is  terminated.  The  kink- 
wave  coordinates,  XWAVEP(S),  YWAVEP(S),  and  ZWAVEP(S),  are  printed  along 
with  the  corresponding  arrestment,  time  TIME  for  both  port  and  starboard 
kink  waves. 


d.  The  remaining  outputs  are  recorded  on  magentic  tape  during 
the  program  calculation  run  and  are  printed  after  the  calculations  are 
completed.  The  maximum  values  are  printed  as  explained  in  Section  III, 
paragraph  E. 

#  e.  The  vehicle  motion  parameter  values  that  are  printed 
along  with  the  arrestment  time  are: 


Speed 

-  VBLNS2 

Runout 

-  VRUN02 

Rolling  friction 

-  VROLL 

Total  drag 

-  VTBRAG 

Acceleration 

-  VACC2 
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Also  printed  out  In  this  output  are: 

Arresting-hook  load  -  VHOOK 

Port  tape  tension  -  DTENP 

Starboard  tape  tension  -  MENS 

The  port  tape  and  tape-reel  values  that  are  printed  out 

Tape  reel  velocity 
Tape  reel  position 
Tape  wrap  radius 
Tape  on  deck 

Arresting-hook-point-to-aheave  distance 

Pendant  elongation 

Total  elongation 

Tape  elongation  factor 

Tape  tension 

Tape  on  reel 

Total  reel  inertia 

Tape-reel  acceleration 


are: 

RSPPR 

RPOP 

ROWP 

TLENP 

DLENP 

PELP 

DELP 

TELFP 

DTENP 

RLTP 

RTINP 

RACPR 
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VII  DISCUSSION 


A.  Run  Verification:  One  of  the  most  important  parameters  associated 
with  arrestment  test  programs  is  arresting-hook  load.  It  is,  therefore, 
used  as  a  basis  of  comparison  to  determine  the  general  accuracy  of  a  com¬ 
puter  mathematical  model  simulation.  During  the  development  stages  of  the 
mathematical  model,  arresting-hook-load  values  obtained  from  a  computer 
run  are  plotted  against  the  arresting-hook-load  values  of  actual  test 
data  to  give  an  indication  of  how  well  a  particular' program  design  logic 
change  has  affected  the  predicted  results.  Other  parameters  are,  of  course, 
analyzed  but  the  arresting-hook  load  serves  as  an  index  of  the  overall  per¬ 
formance  of  the  model. 


B.  Current  Status  of  Program  Performance 


1.  At  the  present  stageiof  development,  the  program  computes  all 
required  tape -reel  and  vehicle  parameters  for  ON- or  OFF-CENTER  arrestments 
from  the  time  of  pendant  pickup  until  the  vehicle  comes  to  a  stop.  How¬ 
ever,  program  logic  at  this  point  of  development  is  not  complete.  The 
basic  logic  pieces  are  present  and  do  account  for  a  complete  parameter 
output,  but  a  need  exists  for  refinement  in  existing  logic  pieces  and  also 
logic  supplementation  to  further  improve  the  accuracy  of  the  model.  The 
problem  that  exisfs  is  the  time  of  peak  value  occurrence  in  arresting- 
hook-load  and  tape-tension  plotted  time  histories.  Peak  values  are  within 
+  5%  of  the  corresponding  actual  peaks;  however,  the  model  calculates  them 
to  occur  at  earlier  times  than  actual  data  indicated  they  occur. 


2.  The  current  logic  piece  that  is  suspected  of  being  the  prime 
cause  for  the  program's  generation  of  peak  values  being  off  in  time  is 
the  determination  of  tape  tension.  Tape  tensions  are  calculated  as  func¬ 
tions  of  purchase-tape  elongation.  An  equation  is  derived  that  expresses 
tension  in  terms  of  percent  elongation  of  the  tape  by  fitting  a  curve 
through  actual  load -percent  elongation  data.  The  problem  is  that  with  nylon 
tape  the  relationship  between  tension  and  percent  elongation  is  hysteretic 
in  nature.  If  a  system  could  be  developed  to  express  the  load-percent 
elongation  relationship  more  adequately,  the  arresting -system  physical 
parameters  could  be  described  more  accurately  by  the  program. 


3,  The  following  is  a  list  of  arrestment  phenomena  that,  up  to 
this  point  in  the  development  of  the  model,  are  not  yet  included: 

a.  Tape -reel  stack  slip 

b.  Angled  arrestments  where  tus  aircraft  or  deadload  does 
not  engage  the  pendant  perpendicularly. 

c.  Arresting-hook ‘slip  on  pendant  during  OFF-CENTER  or 
angled  arrestments 
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f 


d.  Two-block  logic  which  will  provide  for  parameter  values 
after  two-block  condition, 

e.  Effect  that  tape/pendant  connector  has  on  kink -wave 

motion 


f.  Pre-tensioning 

g.  Logic  that  would  account  for  energy  and  momentum  of  the 
arresting  system  with  time 
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VIII  PROGRAM  USE 

A.  Program  Deck  Construction:  The  following  illustrates  the  setup 
of  a  program  card  deck  for  a  normal  run  on  a  CDC  6600  computer: 


NATF-EN-1120 


B.  Input  Data  Preparation:  Program  input  consists  of  six  cards  which 
are  divided  into  two  sections:  arresting-gear  configuration  and  event 
information.  The  cards  are  placed  at  the  end  of  the  card  deck  and  are 
separated  from  the  program  by  a  7,  8,  9  multiple  punched  card.  The  input 
cards  variable  name  location  and  format  are  tabulated  below: 


INPUT 

CARD 


ORDER 

VARIABLE  NAMES 

FORMAT 

ARRESTING-  [ l 

DKSPAN, 

SPLIT?,  SPLITS,  PENLEN,  PENX, 

POFFC 

7F10.4 

GEAR  1  • 

PENDEN 

CONFIGURA-  { 2 

TTHICK, 

TWGHT,  TOLENP.TOLENS,  Cl,  C2, 

C3 

7F10.4 

taoN  1 3 

RINP,  RINS,  BRAKEC,  REUB,  BRAKBX,  TWIDE 

6F10.4 

INFORMATION  1  4 

'PENMOD 

1F12.2 

EVENT  js 

VELENG, 

WEIGH,  VTHRUS,  WINDKT,  BARO, 

TEMPA 

6F10.4 

INFORMATION  \6 

DELTIM, 

VCOAD,  DOFFC,  HOOKHI,  fHKLEN 

5F10.5 

Definitions  of  variable  names  can  be  found  in  Appendix  B.  Input  values 
should  be  obtained  from  the  test  engineer  and/or  the  Computer  Division. 

C.  Control  Card  Preparation:  There  are  five  control  cards  which  pre- 
ceed  the  program  in  the  card  deck  and  are  separated  from  it  by  a  7,  8,  9 
multiple  punched  card. 

1.  Charge  Card:  This  card  is  U3ed  for  accounting  purposes: 

/  $  CHARGE,  aaaaa  C1C2  -  UUU 
/  aaaaa  *  five -digit  charge  number 
CIC2  ■  two  check  digits 

■  indicates  subcharge  number 
UUU  *  three-digit  subcharge  number  (optional) 


2.  Scope  Job  Card:  The  parameters  of  this  card  describe  the  job's 
priority,  time  limits,  memory,  and  peripheral  equipment  requirements.  All 
parameters  are  octal  values: 


Job  Name  (CMfl,  Tt,  10,  Pp) 


Job  Name  =  To  provide  user  identification  of  this  job. 

CMfl  =  Loading  field  length,  which  specifies  the  amount 
of  memory  needed  to  obtain  a  control  point  and 
load  the  program  into  memory. 

Tfc  =  Central  processor  time  limit  for  the  job  in  seconds; 

a  maximum  of  five  octal  digits. 

TO  •=  Inpuc/Output  time  limit  in  seconds. 

Pp  =  Priority  level  (P.)  at  which  job  enters  the  system. 
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Run  Card;  The  FORTRAN  compiler  is  called  by  thiacard: 

f - - — “ - - - 

RUN  (cm,  fl,  bl,  if,  of,  rf,  lc) 

Only  two  parameters  are  required,  cm  and  lc,  the  rest  are 
omitted 

cm  *  Compiler  mode  option 

lc  ■  Line-limit  (octal)  on  the  output  file  of  the  object 
program 


4.  Set  zero  card:  This  card  sets  equal  to  zero  all  storage 
locations  containing  variable  names  and  arrays: 

^/sET  (0) 

5.  Execute  Card:  This  card  loads  the  program  into  memory  and 
begins  execution: 

/  LGO/RFL 

/  LGO  *  Load  and  go 

RFL  *  Request  field  length  -  an  option,  that  allows  optional 
usage  of  central  memory 


For  more  detailed  information  on  preparing  the  control  cards,  see 
references  (c)  and  (d)  or  consult  a  programmer  from  the  Computer  Division, 


V 
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IX  FUTURE  WORK 

A.  Continued  development  on  the  model  will  be  directed  to  the  fol¬ 
lowing  areas: 

a.  Hysteresis  representation  of  the  tape  modulus  of  elasticity 

b.  On-center  and  off-center  angled  arrestments 

c.  Tape  stack  slip 

B.  Also,  at  present,  the  program  outputs  data  in  tabulation  form. 
Plot  capabilities  should  be  added  to  the  program. 
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ARRESTING- GEAR  CONFIGURATION  (INPUT  VALUES) 


PARAMETER  NAME 

FORTRAN 

SYMBOL 

VALUE 

UNITS 

DECK  SPAN 

DKSPAN 

164.000 

FT 

PENDANT  LENGTH 

PENLEN 

154.000 

FT 

PENDANT  MODULUS  OF  ELAS. 

PENMOD 

1872000000.000 

LB/FT- SQ 

PENDANT  AREA 

PENX 

.009 

FT-SQ 

PENDANT  DENSITY 

PENDEN 

17.890 

LB  SEC2/FT4 

PORT  SPLIT  DISTANCE 

SPLITP 

49.000 

FT 

PORT  TAPE  LENGTH 

TOLENP 

840.000 

FT 

PORT  REEL  INERTIA  (METAL) 

RINP 

114.000 

SLUG  FT-SQ 

BRAKE  CONSTANT 

BRAKEC 

12.310 

NONE 

TAPE  LOAD-%  ELONG.  COEFF . 

Cl 

1.709 

NONE 

TAPE  LOAD-%  ELONG. COEFF. 

C2 

85.367 

NONE 

TAPE  LOAD-%  ELONG.  COEFF. 

C3 

4019.790 

NONE 

PENDANT  OFF-CENTER  DISTANCE 

POFFC 

0.000 

FT 

PURCHASE  TAPE  THICKNESS 

TTHICK 

.344 

IN. 

PURCHASE  TAPE  WEIGHT/FOOT 

TWGHT 

1.000 

LB/FT 

fUJICHASE  TAPE  WIDTH 

TWIDE 

8.000 

IN. 

PURCHASE  TAPE  AREA 

TAREA 

2.752 

IN.-SQ 

STBD  SPLIT  DISTANCE 

SPLITS 

49.000 

FT 

STBD  TAPE  LENGTH 

TOLENS 

840.000 

FT 

STBD  REEL  INERTIA  (METAL) 

RINS 

114.000 

SLUG  FT-SQ 

BRAKE-  EXPONENT 

BRAKEX 

2.000 

NONE 

HUB  RADIUS 

RHUB 

9.000 

IN. 

EVENT 

INFORMATION 

(INPUT  VALUES) 

VEHICLE  ENGAGING  SPEED 

VELENG 

150.000 

KN 

VEHICLE  WEIGHT 

VWEIGH 

53000.000 

LB 

VEHICLE  THRUST 

VTHRUS 

0.000 

LB 

HEAD  WIND 

WINDKT 

0.000 

KN 

BAROMETRIC  PRESSURE 

BARO 

30.000 

IN.  HG  A 

AMBIENT  TEMPERATURE 

TEMPA 

68.000 

DEG  F 

CALCULATION  INCREMENT 

DELTIM 

.001 

SEC 

AERO  DRAG  COEFFICIENT 

VCOAD 

.041 

LB  SEC2/FT2 

DECK  OFF-CENTER  DISTANGE 

DOFFC 

0.000 

FT 

ARRESTING-HOOK  HEIGHT 

HOOKHI 

2.313 

FT 

ARRESTING-HOOK  LENGTH 

VHKLEN 

6.167 

FT 

Figure  3  -  Sample  Output  of  Input 
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Figure  4  -  Schematic  of  Runway  Installation  With  Variable  Name  Descriptions 
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Figure  7  -  Arresting-Hook  Load  Determination  (Port  Side  Component  Breakdown) 
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APPENDIX  A  -  MATHEMATICAI  PERFORMANCE  PREDICTION  MODEL  PROGRAM  LISTING 


THIS_PRo8RAM  MAS  REDESIGNED  TO  REFLECT  A  REVISED  APPROACH;  NEW 
LOGIC  INCLUDES  A  3-DIMENSIONAL  KINK  WAVE  FORM, INCREMENTAL  COMPUT¬ 
ATION  OF  THE  KINK  VELOCITY  TO  ACCOUNT  FOR  THE  VARYING  MODULUS  OF 
ELASTICITY  OF  THE  NYLON  T*P£  AND  LOCATION  OF  THE  KINK  (l.E*  IN 
TAPE  OR  STEEL  PENDANT)}  AND  NEW  LOGIC  FOR  KINK2  WHICH  ASSUMES  THE 
FORMATION  OF  A  NEW  KINK  EMINATING  FROM  THE  HOOK  AS  KINKi  IMPACTS 
THE  SHEAVE. 

DIMENSION  DATUMS! 

READ  IN  ALL  DATA  PERTINENT  TO  THE  ARRESTING  GEAR  CONFIGURATION 
JGEOR  A  0 

1  READ  2;  DKSPAN.SPLITP,  SPLITS?  PENLEN)  PENX,  POFFC,  PENDEN 
READ  2;  TTHICK.  TWGHT,  TOLEnP?  TOLENS ,  Cl#  C2,  C3 

READ  3?  RINP#  RINS.  SRAKEC.  RHUS#  BRAKEX,  TWIDE 
READ .4#  PENMOD 

2  FORMAT  (7F10.4) 

3  FORMAT  (6F10.4) 

4  FORMAT  (1F12.J) 

9  FORMAT (sFiOi?) 

READ  IN  ALL  DATA  PERTINENT  TO  EVENT  INFORMATION 
READ  3i  VELENG#  VWEIGMt  VTHRUS.  WtNDKT#  BARO.  TEMPA 
READ. 5?  OELTIN#  VCOAOyOOFFC?  MOOKHI#  VHKLEN 
INITIALIZATION  OF  VARIABLES 
I«L  A  A 

dtenr«dtens*rspprarspsr*o.o 

DELP*DFlS*RACRR«RACSRi  j)  *  o 
VK  INKPAiVK  I  NKS*0 . 0 

VRUNnT4RSpP*RSPS«RACP*RACS*RP0P»RPOS»RELP2»PELS2*RACP2«RAGS2B|;0 
AVHOOK  a  ADTENP  *  ADTENS  AVRCCG  «  A|SPP  *  ARSPS  «ARRCP«ARACSio.O 
XWAVFP  •  YWAVEP  i  ZMKVEP  ■  XWAVES  •  YWAVE8  •  ZWAVEl  ■  Oil 
ICOUNT  a  0.0 
TIME  *  0.0 

VMASS  a  VWEISH  /  32.174 

THDEf  «  DELTIM  /  2;0 

WINDFT  a  WlNDKT  *  1.6178- 

VELNST  a  VELENG  •  1.6878 

VROLt  a  ,00041  •  VELENG  •  VWEIQH 

DRAGS  6  17.3282  •  (  BARO/  (TEMPA  •  499.4  >>  •  VCOAD 

VADRAG  a  (VELNST  NINDFT)M2  •  DRAIN 

VTDRAG  a  VADRAG  ♦  VROLL 

HOOKP  a  HOOKS  *  Ivf 

VHOOK  6  HOQKP  HOOKS 

VACC*VRCC2«<VVHRUS-VM00K*VTDRAG) /VMASS 

VACC0  a  VACC2  /  32vl74 

HYPOTP  a  HYPOTS  «  OtO 

TELFP  a  TELFS  »  8*0 

PELP?  a  OiO 

GAMAYP  a  GAMAVS  *  1*5708 

GAMAXP  a  QAMAIP  a  GAMAXS  a  QAMAZS  «  0.0 

ATTACH  a.  SORT  (VHKLEN  «»2  -  HOOKHI  •  •?> 

DHAP  «  .5  •  DKSPAN  *  DOFFC 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

li 

*2 

13 

14 
19 
14 

17 

18 

1? 

1 1 
n 
*2 

23 

24 

29 
26 
27 
2B 
** 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 


49 

50 
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DH A S  »  ,5  »  DKSPAN  *  ODFFC  5t 

HLENP  *  DHAP  52 

DLENS  k  DHAS  53 

TINP  *  .5  «  (DKSPAN  -  PENLEN  )  ♦  POFFC  54 

TINS  8  ,5  *  (DKSPAN  -  PENLEN  )  -  POTTC  55 

PLENP  *  DHAP  •  TINP  56 

PLFNS  k  DH AS  *  TINS  57 

TLENOP  s  TLENp  ■  SPLITP  ♦  TINP  5B 

TLPNOS  r  TLENS  ■  SPLITS  ♦  TINS  59 

RLTP  «  TOLENP  -  TLENP  AO 

RL-TS  »  TOLENS  -  TLPNS  61 

RNUIP  *  (-TTHICK  -  2»PHU8  ♦  SORT  ((TTHICK  ♦  2*RHUB)#42  ♦  62 

1(  A  •  TTHICK  *  (<  TOLENP  -  TLENP  >  *  !2)/3.14i6  )))7  (2*TTHICK>  63 

ROWP  =  RHUB  ♦(TTHICK  *  RNUlP)  64 

RNU1S  =  (-TTHICK  -  2»RHUB  ♦  SORT  ((TTHICK  ♦  2*RHUB)*#2  ♦  65 

1(  A  *  TTHICK  *  ((  TOLENS  -  TLENS  )  *  12)/3.1416  >>>7  (2»TTHICK>  66 

ROWS  *  RHUB  ♦  (TTHICK  •  RNUIS)  67 

CA  *  TTHIOK  /  6?2832  68 

PEAP  r  PLENP  /  (PENMOO  •  PENX)  69 

PEAS  a  PLENS  /  (PENMOO  •  PENX)  70 

CD  a  (  RHUB  /  12.0  >»«2  71 

CE  a  TWGHT  /  64.348  ‘  72 

LTCONT  =0  73 

TAREA  k  TTHICK  •  TWIDE  74 

10  *  0  75 

Datum ( 32 )  »  o.o  76 

JTIME  k  i  77 

TAPDFN  a  TWGHT/ (  3?-;  17*t  T  AREA/144 . 0  >  >  78 

TEST  FOR  ON  CENTER  ENGAGEMENT  AND  SYMETRICAL  A.G.  CONFIGURATION  79 

ITLGSM  a  0  80 

JSYMVP  3  DHAP*TINp«SPLITP*TOLENP*RINP  8l 

JSYMtfS  .  pHAS4TINS*SPLlTS»T0LENS»RINS  82 

IF(  JSYNVP‘-JSYHVS)19-»18,19  83 

18  IFtGSM  a  i  84 

19  CONTINUE  85 

28  F0RM«T(iHi,///46X,28HARRESTlNG  GEAR  CONFIGURATION)  86 

21  FORMAT  (/28X.7HF0RTRANa52X*7HF0RTRAN)  87 

22  F0RMAT<6X,1AHPARAMETER  NAME 8X  ,  6HS  YM8QL ,  8X ,  5HVALUE ,  3X  ,  5HUNI TS  ,  ilX ,  88 

114WPARKMETER  NAME  ,  8X » 6HSYMB0L  il  8X » 5HVALUE » 3X ,  5HUN I TS7  89 

23  F0RM«T(/10H  DECK  SPAN , j 8X , 6HDKSP AN , IX j F14 , 3 , IX , 2HFTT9X , 24NPENDANT  9q 

10FF  CENTER  DIST.  ,3X»5HP0FFCi'2X,F14.3»lX,2HFT)  *1 

24  FORMAT  (15H  PENDANT  LENGTH , 6HPENLEN , IX , Fi4 . 3 , IX J 2HFT , 9X » 23HPUR  92 

1CHASF  TAPE  THICKNESS. 4X»6HTTHICK.1X»F14. 3, IX, 2HIN)  93 

25  FORMAT  (25H  PENDANT  MODULUS  OF  ELAS.»3X,6HPENM0D,lX»Fl4,3»lX, gNLB/  94 

IFT*S0,-3X,25HPURCHASE  TAPE  WE  I GHT/FOOT  >  2X ,  5HTWGHT  ,  2X  ?F  14 . 3  i  IX ,  5HLB/  95 
2FT>  96 

26  FORMAT  ( 13H  PENDANT  ARFA,l5X»4HDENX>3X»Fl4.3»lX,5HFT-SQ»6X» 19HPURC  97 

1HASE  tape  WIDTH, 8X,5HTWlDEt2)',Fl’.3, IX, 2HIN)  98 

27  FORMAT  (16H  PENDANT  DENSITY '/12X,6HPENDEN,lX,Fl4, 3, llHLB  SEC2/FT  4  •  99 

llX.IRHPURCHASE  TAPE  ARF A , 9X5HT ARE  A , 2X » F14 . 3 , IX , 5H I N*SQ )  100 
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28  FORMAT  ( 20H  PORT  SPLIT  OISTANCE , OX ,6HSPLl TP.1X.F14 , 3 , IX . 2HFT , 9X,  101 


119HST8D  SPLIT  DISTANCE, 8X.6HSPLITS, IX, Fl4, 3, IX, 2HFT)  102 

29  FORMAT  (17H  PORT  TAPE  LENGTH,llX,6HT0LENP,1X,F14.3,iX,2HFT,9X,  103 

116HSTB&  TAPE  L6NGTH,11x,6HT0LENS.1X,F14.3,1X,2HFT)  104 

30  FORMAT  (26«  PORT  REEL  INERTIA  (METAL ) *2X  »4HRInP »3X,Ft4.3»lX,  105 

110HSI‘UB  FT-S0*1X,2?HSTRD  REEL  INERTIA  (METAL) ,2X»4HrINS,3X,  106 


2F14.3,1X,10HSLUQ  FT-SQ)  10? 

31  FORMAT  (15H  BRAKE  CONSTANT, l3x»6HBRAKEC, IX, Fl4.3»iX*4HN0NE,7X»  108 

liAHBRAkE  EXPONENT, 13X,6HBRAKEX, IX, F14, 3, IX, 4HN0NE)  109 

32  FORMAT  (25H  TAPE  LOAD-*ELONG.  COEFF. ,3X,2MCi,5X,F1a;3,iX,4HN0NE,  no 

17X » lOHMUB  RADIUS, 17X,4HRHUB,3X,F14,3»1X,2HIN)  111 

33  FORMAT  (25H  TAPE  LOAD-KELONG.  COEFF. ,3X,2HC2,5X,Fi473.1X,4HN0NE>  112 

34  FORMAT! 25H  TARE  LOAD-(FLONO.  COEFF , , 3X* 2HC3 ,5X»F14 .3, IX, 4MN0NE , 7X ,  113 

125HARRESTMENT  IS  SYMMETRICAL)  114 

35  FORMAT! ///5lX,l7HEVENT  INFORMATION)  115 

36  FORMAT! /32X.14HPARAMETFR  NAME ; lu*  » 6HSYMB0L»8X , 5H VALUE , 3X#5HuNITS )  116 

37  F0RMAT!/29X.25HVEH!CLE  ENGAGING  VELOCITY, 2X* 6HVELENG, lX,Fl4.3iiX,  117 

15HKN0TS) 

38  FORM A T< 29 X *  14 H VEHICLE  WEIGHT,13X.6HVWEIGH,1X.F14.3,1X.3HLBS)  119 

39  F0RMAT<29X,14HVEHICLE  THRUST,13X,6HVTMRUS,1X,F14.3.1X,3HLBS)  120 

40  F0RMAT<?9X,0MHEADNIND,19X,6HWIN0KT,1X<F14.3.1X,5HKN0TS)  1?1 

41  FORMAT! 29X* 19HBAR0METR JC  PRESSURE,8X,4HBAR0,3X,F14.3,lX,7HIN  HC  A)  122 

42  F0RMAT!29X,19HAMBIENT  TEMPERATURE ,  8X  ,5HTEMPA,2X  ,FH73 .  IX  .5HDEG  F>  123 

43  F0RMAT(29X.21HCALCULAT10N  INCREMENT, 6X.6H0ELTIM, IX, Fl4. 3, IX, 3HSEC)  1?4 

44  FORMaT(29x,21HAERO  DRAG  COEFFICIENT, 6X,5HVC0AD,2X,F14. 3, IX, llHl'B  S  125 

1EC2/FT?)  I26 

45  F0RMAT(29X.24HDECK  OFF  CENTER  DISTANCE, 3X,5HD0FFC,2X,F14. 3, 1X.2HFT  127 


46  F0RHATi29X.l5HTAlLH00K  HE  I GHT  U2X ,  6HH00KH  l ,  IX ,  Fl4  i  3  ).1X ,  2HFT ) 

47  FORMAT! 29X,15HTAILH00K  LENGTH, 12X,6HVHKLEN.1X,F14j3;iX.2HFT) 

48  FORMAT ! 25H  TAPE  LOAO-(FLONG.  COEFF  ,  ,3X*2HC3*5X»Fl4.3»lX,4MN0NE*7X, 
129HARRESTMENT  IS  NOT  SYMMETRICAL) 

PRINT  ?0 
PRINT  21 
PRINT  22 

PRINT  23.DKSPAN.P0FFC 
PRINT  24.PENLEN,  TTHICK 
PRINT  25.PFNM0D.THGHT 
PRINT  26, PFNX, TH I OE 
PRINT  27,PENDEN,TAREA 
PRINT  28, SPLITP, SPLITS 
PRINT  29.T0LENP.T0LENS 
PRINT  30.RINS,RINP 
PRINT  3l,8RAKEC,BRAKEX 
PRINT  32,Cl,RHUB 
PRINT  S3.C2 
IF(IFLGSM,EQ.D49,50 

49  PRINT  34, C3 
GO  TO  91 

50  PRINT48.C3 


128 

1?9 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 
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51  PRINT  85  151 

PRINT  86  152 

PRINT  S7^VEL6NG  153 

PRINT  88»VNE!QH  154 

PRINT  39*VTHRUS  155 

PRINT  40«WINDKT  156 

PRINT  U»BARO  157 

PRINT  *2,TEMPA  158 

PRINT  *3,DELTIM  l59 

PRINT  *4,vCOAD  160 

PRINT  JlSiDOFro  161 

PRINT  46>H00KHI  162 

PRINT  *?#VHKLEN  163 

PRINT  Its  164 

88  FORMAT  (  1H1,  32X,  55HPORT  ANO  STARBOARD  RINK  NAVE  COORDINATES  AND  165 

1  VELOCITIES)  166 

PRINT  89  167 

89  FORMAT  (  1  Ho #  llX,  6HP0RT-X?  4X,  6HP0RT-Y,  4X,  6HpORT*Z*  4X,  168 

1  6HVKINKP.  5X »  4HTINE.  15X*  6HSTBD-X.  4X,  6HSTBD-Y.  4X»  6MSTBD-Z.  169 

2  4X,  6WVKINKS,  5X,  4KTIHE)  170 

PRINT  90  ,  XWAVEP,  YNAVEP,  ZNAVEP,  VKINKP,  TIME,  XWAVES,  VMAVES,  l7j 

1  ZWAVES,  VKINKS,  TIME  172 


90  FORMAT  (  10X,  5F10t3,  1 0X*  5F10.3) 

0ATUM(i)«vELNSTSDATUM(9)  >TJM6  $DATUM(i7)»TELFPSDATUM(25>*DLENS 
DATUM(2)»VRUNOTSDATUM(iO)«RSPPR«nATUM(l8)»RLTp  SDATUM(26)«PELS2 
DATUM(3)«vH00M  SDATUM<U)«RpOP  sOATUM(19)»RTInP|DATOM(27)»DELS 
DATUM(4)«VR0LL  SDATUM(12)*R0MP  sDATUM(20)«RACPR1DATUM(28)»TELFS 
DATUM(5)«VTDRAGSDATUMH3)»TLENPSDATUM(21)«RSPSRSDATUM(29)tRLTS 
DATUM(6)»dTENR  |DATUM(14)«0LENPSDATUM(22)«RP0S  SDATOM(30>*RTINS 
DATUM(7)«DTENS  SDATUM(i5)»PELP2SDATUM(23)«R0HS  !DATUM(3l)tRACSR 
DATUM(«)«vACCG  SDATUM<i6)«DELP  SDATUM(24)«TLENS|DATUM(32)*0»0 
WRITE  (7)  DATUM 

MAIN  COMPUTATION  LOOP  BEGINS  AT  STATEMENT  100 
99  CONTINUE 

100  LTCONT  *  LTCONT  ♦  1 

TIME  »  DELTIM  •  LTCONT 

TEST  FOR  PORT  TWO*Bt'OCK 

IF  (  RRUB  .LE,  ROMP  >  I01il008 
TEST  FOR  STBD  TNO40LOCK 

101  IF  (  RMUB  .LE,  RONS  )  102,1008 

TEST  FOR  INSTANTANEOUS  VEHICLE  VELOCITY  »  0 

102  IF  ('  VFLNST  .LE,  OTO  )  1000’,'  103 
STATEMENT  103  BEGINS  PORT  CALCULATIONS 

103  VELNS2  s  VELNST  ♦  THDEL  «  (VACC  ♦  VACC2) 

VRUN02  =  vRUNOT  ♦  THDEL  *  (VELNST  ♦  VELNS2 ) 

RSPP?  s  RSPP  ♦  THDEL  *  (RACP  ♦  RACP2  > 

RPOP?  s  RPOP  ♦  THDEL  *  (RSPP  ♦  RSPP2 ) 

ROWP?  =  ROWP  ♦  CA  •  (RPOP  *  RP0P2) 

RLTP  *  TOLENP  -  TLENP 

TLENP2  =  TLENp  -  ROWP?  *  (RPOP  -  RP0P2  >  /  12.0 
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IF  <  H^POTP/Gt.PLENP)  104,105 

104  JF(T^LPP»LE* .50)  TELFP  *  .50 

ETAPFP  a  -160.635#TELFP**2  ♦  4079800 . 0«TFLFP  -  1919200.0 
CTAPFP  •  SORT  (ETAPEP  /  TApOEN  ) 

TASIOP  ■  ,5  •*  ,667  *  (VELNS2  /  CTAPEP  )  •*  1.33 

VKINKP  a  CTAPEP  *  (SORT  (TaSJGPU  TASIGP)  -(TLENP?  -  TLFNP)  / 

1  OELTIR 
00  TO  106 

105  PEnMOO  «  pEnMoD 

CPENO  .  SORT  (  PENMOO  /  PENQEN) 

PENSJO  a  ,5  *•  .667  •  (  VELNS2  /  CPEN6  )  »4  i,33 

VKINKP  a  CPEND  •  (SORT  (PENSJG)  -  PENSIQ)  -  (TLENPt  -  TLENP  >  / 

1  DELTIIN 

106  XMAVFP  a  XWAVEP  ♦  (VKINKP  •  DELTIM  *  COS  (  GAMAXP  1) 

YNAVPP  a  rWAVEP  -  (VKINKP  •  DELTIM  «  COS  (  GAMAZP)  »COS  (GAMAYP>) 
IMAVPP*  ZWAVEP  *  (VKINKP  *  DELTIM  •  SIN  (  GAMaZP)) 

NK6LPV  a  HOOKHI  -  VHKLFN  •  COS  (1.5708-ATAN  (HOOKHI 7( ATTACH  ♦ 

1  VRU^02 ) ) ) 

HYPOTP  a  SORT  < XWAVF.P*»2»  (VRUN02-YWAVEP)a*2*  (HKELEV-ZWAVEP) 6*2) 
RESUlP  a  SORT  ((DHAP  -  XWAvEp>*»2  ♦  YWAVEP**2  ♦  ZWAVEP«*2) 

DLENP  a  HYPOTP  ♦  RESULP 

DEL*P  a  HYPOTP  ♦  RESUlP  ♦  SPlITP  -  TlENP2  -  PLENP 
IF  (  DELP  .1  f'\  0  >107,108 
18?  DELP  «  0.0 
108  PELP  A  PELP2 

TELFP  a  < DELP  -  PELP  )  /  TLENP2  *  100.0 
IF  (  TfeLFp  .LE.  0  >  109,118 

189  DTENP  a  TELFP  a  PELP2  a  0,0 

118  DTENP  A  Cl  •  TELFP  »*3  ♦  C2  •  TELFP  **2  ♦  C3  •  TELFP 
IF  (DTEnP.LT.  0)  DTENP  a  0 
PELP9  *  PEAP  •  DTENP 

RTINP  4  R I Np  *  CE  *  (TOLENP  -  TLENP2  >*(  CD  ♦  (  ROWP?  /  12.0>**2> 
RACP  a  RACP2 

RACP9  k  (DTENP  *  R0NP2  /  12.0  -  BRAKEC  *  RSPP2»#BRAKEX >  /  RTINP 
HOOKP  k  DTENPaCOS  (  ATAN  (  XWAVEP  /  SORT  ((  HKELEV*ZWAVEP) 

1  «2  ♦  (VRUN02  -YWAVEP)  **2  )>) 

ALL  PORT  CALCULATIONS  DONE  -  TEST  NEXT  FOR  KINK  WAVE  LOCATION 
IF v XWAVEP. GE.DHAP >111, 1112 

111  GAMAXP  a  ATAN  (  SORT  (  HKELEV*»2  ♦  VRUN02*#2)/XWAVEP) 

GAMAYP  a  ATAN  <  XWAVEP  /  VRUN02  > 

GAMA7P  a  ATAN  (  HKELEV  /  SQRT  (  VRUN02»*2  ♦  XWAVEP**2  >) 

XWAVFP  a  0.0 


YWAVFP 

ZWAVFP 


VRUN02 

HKELEV 


HYPOTP  *  0.0 

1112  IFdFLGSM.EQ. 0)112, 1111 
1111  RSPS?  4  RSPP2 
RPOS?  4  RpOP2 
ROWS?  4  ROWP? 

RLTS  *  RLTP 
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TLENS2  a  TLENRf  25l 

VKINkS  «  VK1NKP  252 

XUAVpS  >  XNAVfP  253 

YWAVPS  a  YWAVCP  254 

ZWAVPS  •  ZWAVIP  255 

DLENS  a  DLENP  256 

DELS  ■  DELP  257 

TELFS  •  TELFP  258 

DTENS  «  DTENP  259 

PELS?  a  PELP2  26o 

RACS?  *  RaCP2  261 

HOOKS  •  HOOKP  262 

QAMAXS  a  GAMAXP  263 

GAMAYS  •  GAMAYP  264 

8AHA7S  •  QAHA2P  2*5 

GO  TO  130  266 

C  STATEMENT  112  BEGINS  STARSOaRO  CALCULATIONS  267 

112  RSPS?  •  RSPS  ♦  THOEL  •  (  RACS  ♦  RACS2  )  268 

RPOSP  *  RPOS  ♦  THDEL  ••  (  RSPS  ♦  RSPS2  )  269 

ROWS?  A  ROUS  ♦  CA  6  (  RPoS  «  RPOS2  )  270 

RLTS  ■  TOLENS  -  TLENS  27l 

TLENS2  a  TLENS  -  ROWS2  •  (  RPOS  -  RP0S2  )  /  12.0  272 

IF  (WYPOTS  .GE.  PLENS)  113*114  273 

113  IF(TPLFS.LE..50)TELFS  •  .50  274 

ETAPFS  >  -160.635*TELFS**2  ♦  407*800 .8<*TELFS  -  1919108.0  275 

CTAPFS  a  SORT  (ETAPES  /  TApDEN)  276 

TASIGS  a  .5  *4  ,667  •  (VELNS2  /CTAPES  )  ••  1,33  277 

VKINKS  ■  CTAPES  •(SORT  (TASIGS)-  TASIGS)  -  (TLENS2  -  TLENS)  /  278 

1DELTIH  279 

GO  TO  li5  280 

114  PENMDD  a  PENMOD  2§l 

CPEND  a  SORT  (  PENMOD  /  PENBEN  )  2»2 

PENSIG  a  ,5  ••  *667  *  (  VELNS2  /  CPEND  )  •  *  1,33  2«3 

VKINKS  a  cPENO  •  (SORT  (PENSIG)  -  PENSIG)  -  (TLENS2  ♦  TLENS  *  /  284 

1DELTTM  2«5 

115  XHAVPS  a  XWAVES  ♦  (VKINKS  •  DEL T ) M  *  COS  (  GANAXS  )T  286 

YHAVFS  a  YWAVES  -  (VKINKS  *  DELTlM  #  COS  (GAMAZS)  *  COS  (GAMAYS))  287 

ZWAVPS  a  ZWAVES  -  (VKINKS  •  DELTlM  •  SIN  (  GaMA7S>7  200 

HYPOTS  a  SORT  (XWAVES*»2*(VRUN0?-YWAVFS>**2  ♦(HKeLEV-ZWAVES)»42)  209 

RESULS  i  SORT  ((  DMAS  -  XkAVES>**2  ♦  YWAVES*a2  ♦  ZMAVES**2)  2*0 

DLENS  a  HYPOTS  ♦  RESULS  29l 

DELS  a  HYPOTS  ♦  RESULS  ♦  SPLITS  -  TLENS2  -  PLENS  292 

IFID^LS  »LTi  8  )  116.117  293 

It*  DELS  a  0.0  294 

i!7  PELS  a  pELSJ  295 

TELFS  a  (DELS  •*  PELS  )  /  TLENS2  •  106.0  296 

IF  (TSCFS  .LEi  0  )  118; 119  297 

118  DTENR  A  TELFS  a  PEL'82  a  0,0  298 

119  DTENR  *  Cl  *  TELFS  1*3  ♦  C2  *  TELFS  **2  ♦  C3  *  TELFS  299 

IF  (DTInS-.LT.  0)  DTENS  a  0  300 
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120 


ISO 


131 

132 


26$ 

261 

262 

263 

264 

265 
26$ 
287 
208 


PELS?  6  PEAS  4  DTENS 

RUNS  *  RINS  *  CE  #  (  TOKENS  -  TLENS2>*<  OD  ♦  (R0NS2  T  12?0  >#*2) 
RACS  a  RACSJ 

RACS?  A  (  DTENS  •  RONS?  /  12.8  -  BRAKSC  •  RSPS2**  BRAKEX )  /RTINS 
H06KR  A  DTENSiCOS  <  ATAN  (  XNAVES  /  SORT  ( (HKELEV-jNAVES) 

1  »«2  ♦  (VRUNOa-YWAVESf  »«2  )>) 

ALL  STfiD  CALCULATIONS  DONE  /  TEST  NEXT  FOR  KINK  NAVE  LOCATION 
IF  (  XNAVES  .OE.  ONAS  )  1207130 

QAHAXS  a  ATAN  (  SORT  (  HKELEV*«2  ♦  VRUN02*«2>/XNAVES) 

OAMAYS  *  ATAN  (  XNAVES  /  VRUN02  > 

0AHA7S  a  ATAN  (  MKEL'EV  t  SORT  <  VRUNO?**?  ♦  XHAVES*. 2  >) 

XHAVf-S  a  0.0 

YMAVPS  *  VRUN02 

ZNAVfcs  •  HKELEV 

MYPOTS  *0.0 

STATEMENT  430  BEGINS  FINAL  CALCULATIONS 

VHOOK  «  HOOKP  ♦  NOOKS 

VFRIRT  •  .0002666  •  VELNS2 

IF  (.VPRICT  .L*E.  .018  )  1317132 

VFRIfJT  •  .018 

VROLI;  a  VFRICT  •  VNSIGH 

VABRA8  a  (  VELNS2  *  NfNDFT)V*2  •  DRAGK 

VTDRAG  a  yADRAG  ♦  VROLL 

VACCG  •  VACC3  /  327174 

VACC  a  VACC3 

VACC?  a  (  VTNRUS  -  VHOOK  -  VTDRAG)/  VMASS 

RP8P  a  RPOPa 

RPOS  a  RPOSJ 

RONP  a  RONPa 

RONS,  a  RONSa 

TLENP  a  TLENP2 

TLENS  a  TLENS2 

VEUNST  a  VELNS2 

VRUNOT  a  VRUN02 

RSPP  a  RSPPa 


RSPS.  a  RSPSa 
RSPPR  a  RSPP  • 
RSPSR  a  RSPS  • 
RACPR  a  RACP  * 
RACSR  a  RACS  a 


9.5493 
9.5493 
572»9578 
57279578 

SAVE  MAXIMUM  VALUES  BEGINS  AT  STATEMENT  200 
IF  (.AVHOOK  iLT.  VW0OM  2015*202 


AVNOOK  a  VHOOK 
IF  (  ABtENP  .LIT. 
ADTENP  a  DTENR 
IF  <  AOTENS  .LT, 
ADTENS  a  dTENS 
IT  <_AVaCCG  .QT. 
AVACfiG  a  VACCG 
IF  (  ARSPP  -LT. 


DTENPi)  2035*204 
DTENS)  2O5?206 
VACCG)  2075*208 
RSPPR)  209{*210 


S  ATHOOK 
1  ATDENP 
S  ATDENS 
S  ATVACG 


,  ENDS  AT  215 
«  TIME 
a  TIME 
a  TIME 
a  TIME 


301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 
3)5 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 
33$ 

331 

332 

333 

334 

335 

336 

337 

338 

339 
34$ 

341 

342 

343 

344 

345 

346 

347 

348 

349 
35$ 
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269  ARSPP  a  RSPPR  S  ATRSPP  a  TIME  S5l 

2ie  rr  (.aiIsps  .ke.  rspsri  211,212  352 

211  ARSP*  «  RSPSR  S  ATRSPS  a  TIME  353 

212  rr  <_A(*ACP  . kT «  RACPR)  213,214  354 

213  ARACP  «  RACPR  *  ATRACP  a  TIME  353 

214  If  (  . AflACS  .UT.  RACSR)  2l5,2l6  356 

215  ARACS  •  RACSR  S  ATRACS  a  TIME  357 

C  MAIN  UftOP  OUTPUT  AND  TAPE  STORE  LOGIC  BEGINS  HERE  358 

216  CONTINUE  359 

IFtTTMf .EQ.OIGO  TO  300  360 

IGL  «  tGL  4  1  361 

IF(IGL:eQ. 10)402, 400  362 

402  IGL  *  6  363 

300  DATUM(l)avELNSTSDATWM(9)  aTlME  jDATUM(l7)aTELFP$DATUM(25)«DLENS  364 

DATUM<2)«VRUNOTSDATUM(iO)aRSPPRSDATUM(18)aRLTP  *DATUM(26)apELS2  365 

DATUM < 3 ) « vHOOK  IOATUM(ll).RpOp  SDATUM{19)iRTlNPSDATWM<27)«DELS  366 

DATUM ( 4 ) ayROLL  1DATUMI j2)aR0WP  SOATUM(20 )aRACPRSDATUM(28)«TELFS  367 

DATUM<9)svTDRAGSDATUM<ij3)aTL6NPjnATUM<21)aRSPSRSDATUM<29>aRLT$  360 

DATUM ( 6 ) =DTENP  SDATUM(i4)a0LENPinATUM(22)aRP0S  SDATBM(30)*RTINS  369 

DATUM ( 7 ) «dTENS  SDATUM<i5)aPELP2SDATUM<23)aR0WS  SDATUM(31)*RACSR  370 

DATUM ( fl ) ivACCG  SDATUM(lO)aDELP  *DATUM(24)aTLENSSDATHM(32)40.0  3?1 

WRIT?  <7)  DATUM  372 

400  IF(TTME.LE.3.0)301,304  373 

301  IF(LTCbNT-10aJTIMET304i302»302  374 

302  PRINT  90  ,  XWAVEP,  YWAVEP,  ZWAVEP,  VKINKP,  TIME,  XWAVES,  YWAVES,  3?3 

1  ZWAVlES,  VKINKS,  TIME  370 

303  JTIMF  i  JTlME  ♦  1  377 

304  GO  TO  99  3?8 

C  MAIN  COMPUTATION  LOOP  DONE  FINAL  OUTPBT  BEGINS  WITH  STATEMENT  1000  379 

1000  REWIND  7  3«0 

c  print  Out  maximum  values  of  parameters  and  time  of  occurance  381 

PRINT  1001  382 

1001  FORMAT  (43X.7HMAXIMUM,  24X>  4HTIME  >  3*3 

PRINT  1002*  AVH00KV  ATHOOK  384 

1002  FORMAT  (10H  HOOKLOAD*  2*Xi  F12.4*  3H  LB,  13X.  Fl0v4J  4H  SEC)  383 

PRINT  1003*  ADTENP;  ATDENP  380 

1003  FORMAT  (19H  TAPE  TENSION'iPORT.lPX^lZ.A.SH  LB,13X,Pl0.4,4H  SEC)  387 

PRINT  1004.  ADTENS,.  AfDENS  388 

1004  FORMAT  (19H  TAPE  TENSI0NiSTBB,l9X,Fl2.4,3H  LB.13X,Pl*.4,4H  SEC)  389 

PRINT  1005.  AVACCG,  ATVACG  3’0 

1005  FORMAT  (22H  VEHICLE  DECELERATION,  16X,  F12.4,  4H  G  S,  12X,  FlO.4,  391 

1  4H  SEC)  392 

PRINT  1006.  ARSPP ,  ATRSPP  393 

1006  FORMAT  (17H  REEL  SPEED-PORT, 21X , F12 . 4 , 4H  RPM,12X *FiO ,4 ,  4H  SEC)  394 

PRINT  1007.  ARSPS,  ATRSPS  395 

1007  FORMAT  (17H  REEL  SPEED-ST8D, 21X , F12 • 4 » 4H  RPM,12X*Fl0.4,  4H  SEC)  390 

PRINT  1008,  ARACP,  ATRACP  397 

1008  FORMAT  (24H.  REEL  ACCELERATION-PORT,  14X,  F12.4,  7H  RPM  SO,  9X?  398 

1  FlO  » 4  ,  4H  SEC)  399 

PRINT  1009.  ARACS,  ATRACS  400 
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1109  FORMAT  (24H  REEL  ACCELERaTION-STRD*  J4X*  Fl2.4>  7H  ROM  90,  9Xf 
1  FiO*4?  4H  SEC) 

1110  JTTMfc  i  1 
1056  JCOUNT  *  0 

PRINT  jf 052 

1058  FORMAT  ( 9H1  VEHICLE* 

PRINT  1094 

1054  FORMAT  (106H0  VELOCITY  RUNOUT  HOOKLOaO  ROLL*  FRJCt  fOT 
1AL  DRAG  TENSION/P  TENSION/S  ACCELERATION  TIME) 

PRINT  i056 

1056  FORMAT  (106H  FT/SEC  FT  LB  LB 

1  LB  LB  LB  G  SEC, 7/> 

1058  READ  < h  DATUM 

i860  LTCONT  «  pATUM<9)*1000~0 

IF(LTCbNT-10*JTIMET105B >1064*5064 

1064  PRINT  <eo»6*  (  DATUM  (!)>!■  1*9  > 

IF (DATUM! 9 ) J  GE . 2 . 0)1063*1065 

1063  JTTMi?-*  JTIME  ♦  10 
GO  TO  Il068 

1065  JTIME  A  JTIME  ♦  1 

1066  FORMAT  <7F12>2,Fl274»Fl2i 8) 
i860  .TIM  4  'tlMEu*0l 

rF(DATUM(9);GI.TIMTI072*1070 
1070  JCOUNT  •  JCOUNT  ♦  1 

IF  (  JOOUNT  .OE.  50  )  1050(1050 
1072  REMIND  7 
1673  mMM  1 
1074  JCOUNT  *  0 
PRINT  |o76 

1076  FORMAT  (HH1  RORT  SIDE) 

PRINT  5o7B 

1078  FORMAT  (loBH  REEL  REEL  RADIUS  TAPE  HOOK  PT  PEND 

1  TOTAL  ELONO  TAPE  TAPE  ON  TOTAL  REEL) 

1080  F0RMAT1(U6H  SPEED  POSITION  NRAP  ON  DECK  TO  SHEAVE.  EbONG 

1  ECONO  FACTOR  TENSION  REEL  INERTIA  ACC  TIME) 

PRINT  *082  t-  „  -- 

10O2'  format  di6H  rpm  rad  in  ft  ft  ft 

1  pt  %  LB  FT  SLOG-FT  SQ  RPM  SQ  SEC) 

1084  READ. t7)  DATUM 

1086  LTCONT  s  dATUM<9>*1000$0 
IFILTC0NT  w  10*JT IME ) 1084 il090* 1090 

1090  PRINT  lo92»  <  DATUM! I ) j  W10.67  )»  DAfUM(6>,  DATUM(i8)» 

1  DATUM (i9) .  DATUM(20),  DATUH<9> 

IF(D*T0M(9)iG6.2.0)1087<l09l 

1087  JTIM^  4  JTIME  ♦  10 
GO  TO  |o 94 

1192  FORMAT^ <F9»3%i0^3*F9; 3*P872*P9 .3 1 2(F8»3)  *'F8.4(F10«  1  »F8.‘2»F12. 3. 
lFil,3/F8.3> 


A-9 


NATF-EN-1120 


1094  TIM  ?  *TlME-.Ot 

IF<D4Tt)M<9)  j S6. TlMlf  1090 #1096 
1096  JCOUNT  s  JCOUNT  ♦  1 

IF  <  JGOUNT  . QE,  50  >  107411004 
1098  REMIND  7 

JQEOR  i  JGEOR  ♦  1 
2000  ST0P 
END 


451 

452 

453 
$54 
455 
450 
457 
450 
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APPENDIX  B  -  VARIABLE  NAMES,  TYPES,  DEFINITIONS,  AND  UNITS  (NAMES  ENDING  WITH  THE  LETTER  "P"  REPRESENT  PORT 
SIDE  VALUES;  THOSE  ENDING  IN  "S"  REPRESENT  THE  CORRESPONDING  STARBOARD  VALUES) 


Type 

Variable 

Name 

In- 

PUt 

Out- 

put 

Inter¬ 

me¬ 

diate 

Definition 

Units 

ADTENP(S) 

X 

Maximum  tape  tanalon 

Lb 

ARACP(S) 

X 

Maximum  raat  angular  accalaratlon 

Rad/Sec2 

ARSPP(S) 

X 

Maximum  real  angular  velocity 

Rad/Sec 

ATDENP(S) 

X 

Tima  at  which  AVHOOK  occurs 

Sac 

ATHOOK 

X 

Tima  at  which  ADTENP(S)  occurs 

Sec 

ATRACP(S) 

X 

Tima  at  which  ARACP(S)  occurs 

Sac 

ATRSPP (S) 

X 

Tima  at  which  ARSPP(S)  occurs 

Sec 

ATVACG 

X 

Time  at  which  AVACCG  occurs 

Sec 

AVACCG 

X 

Maximum  vehicle  acceleration 

G 

AVHOOK 

X 

Maximum  vehicle  arresting-hook  axial  load 

Lb 

BARO 

X 

Barometric  pressure 

In.  Hg 

BRAKEC 

X 

Tape  reel  brake  constant 

None 

BRAKEX 

X 

Water  brake  exponent 

None 

Cl 

X 

Coefficient  for  tape  tension  equation 

None 

C2 

X 

Coefficient  for  tape  tension  equation 

None 

C3 

X 

Coefficient  for  tape  tension  equation 

None 

CA 

X 

Tape  thickness 

2it 

In. 

CD 

X 

^Hub  radius^2 

Ft2 

CE 

X 

Tape  welaht  per  foot 

2g 

Lb-Sec2 

Ft 

CPEND 

X 

Stress  propagation  velocity- pendant 

Ft/Sec 

CTAPEP(S) 

X 

Stress  propagation  velocity-tape 

Ft/Sec 

DATUM(32) 

X 

Linear  array  dimension 

None 

DELP(S) 

X 

Total  elongation  (tape  and  pendant) 

Ft 

DELTIM 

X 

Incremental  time  for  each  calculation 

Sec 

DHAP(S) 

X 

Distance  -  sheave  to  point  of  engagement 

Ft 

DKSPAN 

X 

Distance  between  runway-edge  sheaves 

Ft 

DLENP(S) 

X 

3D  distance,  hook-point-to-klnk-to-sheave 

Ft 

DOFFC 

X 

OFF-CENTER  engaging  distance  (starboard  +) 

Ft 

DRAGK 

X 

Drag  factor 

Lb-Sec2 

Ft2 

DTEHP(S) 

X 

Tension  In  tape  and  pendant 

Lb 
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Variable 

Name 

-  Type 

Definition 

Units 

In- 

£U£_ 

Out- 

£ut_ 

tnter- 

dlats 

EIAPEP(S) 

X 

Tape  modulus  of  elasticity 

Lb/Ft2 

GAMAXP(S) 

X 

Angle,  Initial  pendant  posltlon-to-sheave-to-klnk 

Deg 

GAMAYP(S) 

X 

Angle,  Intersection  of  line  of  engageaent  and  projection  of 

Deg 

line  from  aheave  to  kink 

GAMAZP(S) 

X 

Angle,  klnk-to-sheave-to-projectlon  of  kink  on  deck 

HKELEV 

X 

Hook  point  elevation  at  any  Instant 

Ft 

HOOKHI 

X 

Elevation  of  hook  attach  point  on  vehicle 

Ft 

HOOKP(S) 

X 

Component  of  arreatlng-hook  axial  load 

Lb 

HYPOTP(S) 

X 

3D  distance,  hook  point  to  kink 

Ft 

ICOUNT 

X 

Output  printer  line  counter 

None 

IFLGSM 

X 

Engagement  symmetry  flag 

None 

IGL 

X 

Output  generation  Increment  counter 

None 

JCOUNT 

X 

Line  counter 

None 

JGEOR 

X 

Program  sequential  run  Indicator 

None 

LTCOh’T 

X 

Calculation  Increment  counter 

None 

PEAP(S) 

X 

Used  to  segment  pendant  elongation  equation 

1 

Lb-Ft3 

PELP(S) 

X 

Pendant  elongation 

Ft 

PENDEN 

X 

Pendant  density 

Lb-Sec2 

Ft* 

PENLEN 

X 

Total  unelongated  pendant  length  (eye-to-eye) 

Ft 

PENMOD 

X 

Modulus  of  elasticity  of  the  pendant 

Lb/Ft2 

PENSIG 

X 

Transverse  Impact  equation  approximation  for  pendant 

None 

PENX 

X 

Pendant  cross-sectional  area 

In.2 

PLENP(S) 

X 

Pendant  length  (unelongated) 

Ft 

POFFC 

X 

Pendant  off-center  distance  (starboard  +) 

Ft 

RACP(S) 

X 

Reel  angular  acceleration 

Rad/ Sec2 

RACP(S)R 

X 

Reel  angular  acceleration 

Rev/Mln2 

RESULP(S) 

X 

3D  distance,  kink  to  sheave 

Ft 

RHUB 

X 

Tape  reel  hub  radius 

In. 

RIKP(S) 

X 

Tape-reel  Inertia  (metal) 

Slug-Pt2 

RLTP(S) 

X 

Length  of  tape  on  reel  (not  on  deck) 

Ft 

mie(s) 

X 

Initial  number  of  tape  wraps  on  reel 

None 
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_ Type 

Inter- 

Variable  In-  Out-  me- 

Nama  put  put  dlate  Definition  Pnlte 


ROWP(S) 

X 

Tape  outer  wrap  radius 

In. 

RPOP(S) 

X 

Reel  angular  position  (beginning  of  time  increment) 

Rad 

RPOP2(S2) 

X 

Reel  position  (end  of  tine  Incranent) 

Rad 

RSPP(S) 

X 

Reel  angular  velocity 

Rad/Sec 

RSPP(S)R 

X 

Real  angular  velocity 

Rev/Min 

RTINP(S) 

X 

Total  reel  Inertia 

Slug-Ft2 

SPLITP(S) 

X 

Split  distance  (fr<®  runway  sheave  to  tape  real) 

Ft 

TAPDEN 

X 

Tape  density 

TAREA 

X 

Tape  cross-sectional  area 

In.2 

TASIGP(S) 

X 

Transverse  impact  equation  approximation  for  tape 

None 

TELFP(S) 

X 

Tape  elongation  factor 

X 

TEMPA 

X 

Temperature  (ambient) 

°  F 

THDEL 

X 

mxw 

2 

Sec 

THETAP(S) 

X 

Kink  angle 

Deg 

TIME 

X 

Absolute  accumulative  time 

Sec 

TINP(S) 

X 

Initial  length  connector  to  sheave 

Ft 

TLENOP(S) 

X 

Initial  tape  on  deck 

Ft 

TLENP(S) 

X 

Tape  on  deck  at  beginning  of  time  Increment 

Ft 

TLENP2(S2) 

X 

Tape  on  deck  at  end  of  time  increment 

Ft 

TOLENP(S) 

X 

Total  purchase-tape  length 

Ft 

TTHICK 

X 

Purchase-tape  thickness 

In. 

TOGHT 

X 

Purchase-tape  unit  weight 

Lb/Ft 

•WIDE 

X 

Purchase-tape  width 

In. 

VACG 

X 

Vehicle  acceleration  at  the  beginning  of  the  time  increment 

Ft/Sec2 

VACC2 

X 

Vehicle  acceleration  at  the  end  of  the  time  increment 

Ft/Sec 

VACCG 

VADRAG 

X 

X 

Vehicle  acceleration 

32.17 

Vehicle  aerodynamic  drag 

G 

Lb 

VCOAD 

X 

Vehicle  coefficient  of  aerodynamic  drag 

None 

VELENG 

X 

Vehicle  engaging  speed 

Kn 

VELNST 

X 

Vehicle  speed  at  beginning  of  time  increment 

Ft/Sec 

VELNS2 

X 

Vehicle  speed  at  end  of  time  increment 

Ft/Sec 

VFRICT 

X 

Coefficient  of  rolling  resistance 

None 

VHKLEN 

X 

Vehicle  arresting-hook  length 

Ft 

VHOOK 

X 

Vehicle  arresting-hook  axial  load 

Lb 
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VKINKP(S) 


Inter- 

In-  Out-  me- 
put  put  dlate 


Kink-vave  velocity 
Vehicle  meee 

Vehicle  rolling  reeletence  load 

Vehicle  runout  at  beginning  of  tine  lncremont 

Vehicle  runout  at  the  end  of  tine  Increment 

Vehicle  total  drag 

Vehicle  thruat 

Vehicle  weight 


X  Head-wind  velocity 
Head-wind  velocity 


XWAVEP(S) 


X  Kink  X-dlatance  from  point  of  engagement 


YWAVEP(S) 


X  Kink  Y-dietence  from  point  of  engagement 


ZWAVEP(S) 


X  Kink  Z-dlatance  from  point  of  engagement 


APPENDIX  C  -  GENERAL  FLOW  CHART 


4 
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APPENDIX  D  -  DETERMINATION  OF  KINK-WAVE  VELOCITY 

1.  The  location  of  the  kink  wave  with  time  is  an  important  factor  in 
the  mathematical  model  design.  The  kink  wave's  position  on  the  runway 
determines  the  overall  pendant/tape  geometry  which  governs  the  operation 
of  the  program.  The  change  in  runway  geometry  directly  affects  the  two 
most  important  arrestment  parameters:  tape  tension  and  arresting-hook 
load. 

2.  The  location'  of  the  kink  wave  is  determined  by  its  velocity  alone 
because  its  path  is  predetermined  by  model  design.  In  this  model,  the 
kink  wave,  upon  hook  impact  with  the  pendanc,  follows  a  straight  line 

path  from  the  position  of  hook  impact  to  the  runway-edge  sheave  (kink  one) . 
Upon  impact  on  the  sheave,  it  then  follows  a  new  path  which  is  a  straight 
line  from  the  position  of  the  hook  at  sheave  impact  to  the  runway-edge 
sheave  (kink  two) .  The  establishment  of  kink-wave  paths  is  illustrated 
below: 


3.  The  approach  to  the  problem  of  determining  kink-wave  velocity  is 
based  on  information  presented  in  "Cable  Dynamics"  by  F.  0.  Ringleb 
(reference  (a)  of  basic  text),  and  the  method  of  solution  presented  here 
follows  from  his  work. 

4.  It  is  assumed  that  the  arresting  hook  of  the  aircraft  engages  the 
pendant  perpendicularly  and  results  in  transverse  impact.  A  triangle¬ 
shaped  deformation  is  formed  in  the  pendant  due  to  impact.  A  stress 
wave  propagates  along  the  pendant  toward  the  sheaves  with  a  finite 
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velocity  ahead  of  and  faster  than  the  deformation.  The  velocity  at  which 
the  kink  wave  propagates  along  the  pendant  while  this  is  happening  is  es¬ 
tablished  below? 

NOTE:  Ringleb's  analysis  applies  only  to  an  infinitely  long  cable.  The 
assumptions  made  in  the  model  design  are;  19  The  calculation  of 
the  kink-wave  velocity  in  the  purchase-tape  media  follows  the 
same  principles  that  pendant  calculations  are  based  on  and  2,  the 
equations  of  motion  described  are  used  even  after  kink  one  reflects 
off  of  the  runway-edge  sheave  and  the  pendant/ tape  configuration 
no  longer  appears  to  be  like  an  infinitely  long  cable. 

Letting: 


E  *  modulus  of  elasticity  of  pendant  (tape) 
p  *  mass  density  of  pendant  (tape) 
a  *  stress 
ao  *  initial  stress 

Then  from  the  classical  theory  of  a  vibrating  string,  the  longitudinal 
wave  velocity  C  is: 


/E* 

°-/7 

and  the  transverse  wave  velocity  C  is 

G  -  /f  . 

G  is  the  velocity  of  stress  propagation  in  the  medium  of  concern,  and  C 
is  the  velocity  that  the  deformation  propagates  with  respect  to  a  par¬ 
ticular  defined  mass  point  in  the  media  (pendant)  at  time  of  impact 
(t  *  0).  The  following  diagram  shows  a  segment  of  the  pendant  impact: 


D-2 


NATF-EN-1120 


In  the  preceding  diagram,  the  motion  of  a  point  mass  TJ  of  the  pendant  at 
time  0  is  analyzed  to  define  kink-wave  motion.  After  impact,  t  >0,  this 
point  mass  is  designated  Q  and  marks  the  location  of  the  kink  wave  on  its 
way  toward  the  sheave.  R  represents  the  point  that  the  stress  has  reached 
at  t  >  0._  The  point  Q  moves  toward  the  right  with  the  transverse  wave 
velocity  C  while  the  pendant  segment  QR  moves  toward  the  left  with  a 
velocity  ^  (from  the  longitudinal  impact  formula  a  -  ao/E  =  (i/C).  Thus 
the  velocity  of  the  kink  wave,  W,  is  given  by  the  relationship  W  *  G  -  jj, . 
Now,  C  can  be  expressed  in  terms  of  E  and  a  since 


C 


it  follows  that 


_1_ 

/P 


_C_ 

A 


C  1 

so  that  substituting  for  -±-  into  the  formula  for  £  gives 

A  A 


the  kink-wave  velocity  can  therefore  be  expressed  as 


5.  An  expression  for  the  term  ?  which  represents  an  approxima¬ 
tion  of  the  transverse  impact  formula  (reference  (a)  of  main  text),  is 
given  as, 

g-PQ  /  1  \2/3  /  Vo  \A/3 

E  ”  V  2  /  (  C  / 

where  Vo  is  the  impact  velocity.  So  that  finally,  the  kink-wave  velocity 
for  an  infinitely  long  cable  is: 

NOTE:  The  kink-wave  velocity  in  this  model  is  computed  with  this  equa¬ 
tion  modified  with  an  additive  factor  accounting  for  the  tape 
"feed  in"  velocity  as  the  tape  reel  responds  to  the  acceleration 
forces  during  an  actual  arrestment. 


D-3 


NATF-EN-1120 


APPENDIX  E  -  DETERMINATION  OF  THE  MODULUS  OF  ELASTICITY 

1.  The  value  of  modulus  of  elasticity  for  both  steel  (pendant)  and 
nylon  (tape)  is  required  in  the  program  to  account  for  kink-wave  motion. 
The  value  for  the  pendant  is  a  constant  value  and  is  input  into  the  pro¬ 
gram.  The  value  for  nylon  tape  varies  with  tape  stretch  and  an  expres¬ 
sion  must  be  derived  to  represent  it.  Therefore,  a  mention  of  elasticity 
theory  and  the  method  of  derivation  of  an  appropriate  modulus  equation 
for  nylon  tape  follows. 

2.  The  modulus  of  elasticity  is  a  measure  of  the  effect  that  tension 
has  on  a  material.  Many  materials  exhibit  a  stretching  (deformation)  in 
direct  proportion  to  the  amount  of  tension  applied  (loading) .  This  re¬ 
lationship  is  given  by  Hooke’s  law: 


where 

and 

If 


and 
then 
so  that 


E  »  modulus  of  elasticity  lb/ in.  , 
S  =  stress  lb/in.^, 
e  =  strain  in. /in. 

L  *  original  length  of  material, 

A  =  original  area  of  material, 

P  =  applied  tension  load, 

5  ~  elongation  of  material, 

S  =  P/A  and  e  =  6/1 
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Plotting  P/A  against  6/L,  we  obtain  the  typical  stress-strain  curve: 


Figure  El 

Line  OA  is  linear  and  represents  the  elastic  range  of  a  material.  In 
this  range,  no  permanent  deformations  are  produced  by  stresses.  Hooke's 
law  applies  only  in  this  region.  Curve  AB  represents  the  plastic  range 
of  a  material.  Here,  permanent  deformations  or  sets  occur  upon  applied 
loading.  Point  B  shows  where  the  test  material  fails.  The  slope  of  the 
line  0&  is  the  modulus  of  elasticity  of  the  material.  Stresses  and 
strains  along  this  line  follow  the  "linear  theory  of  elasticity". 

3.  Modulus  of  Elasticity  of  Steel  Pendant:  The  stress-strain  dia¬ 
gram  for  a  steel  pendant  closely  resembles  that  of  Figure  El.  The  value 
of  modulus  used  in  the  program  is  taken  from  data  for  steel  wire  cables 
with  hemp  core  containing  6  strands  with  19  wires  each.  The  value  is 
E  =  13,000,000  lb/ in. ^  (page  13  of  reference  (a)  of  basic  text).  This 
value  varies  considerably  from  that  of  steel--E  *  30,000,000  lb/in. 2-- 
indicating  that  a  pendant  of  this  type  is  less  stiff  than  a  correspond¬ 
ing  solid  element. 
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4.  Modulus  of  Elasticity  of  Nylon  Tape;  Due  to  the  inhomogeneity 
of  the  nylon  threads  and  the  fabric  nature  of  nylon  tapes,  the  stress- 
strain  diagram  representing  the  tensile  test  of  a  piece  of  nylon  tape 
differs  considerably  from  that  of  Figure  El: 


Figure  E2 

Nylon  tape  used  on  the  arresting  gears  are  actually  plastic  in  nature. 

When  a  load  is  applied  to  a  specimen  in  stepped  increments  and  then  un¬ 
loaded,  the  curve  representing  this  action  does  not  return  to  where  it 
originated  but  shows  evidence  of  permanent  tape  elongation.  Upon  re¬ 
peated  loadings  and  unloadings,  loops  are  formed  representing  the  tape's 
inelasticity.  After  the  initial  pull,  the  loops  that  are  generated 
originate  in  approximately  the  same  area.  Also,  the  plastic  deformations 
are  somewhat  predictable  and  therefore  can  be  thoughc  of  as  uniform  plas¬ 
tic  deformations.  For  these  reasons,  a  value  for  modulus  of  elasticity 
of  nylon  tape  can  be  obtained  by  considering  the  uniform  plastic  deforma¬ 
tion  to  be  elastic  in  nature  and  thereby  obtain  an  approximate  value  for 
the  modulus  by  applying  elastic  theory  methods  (modulus  is  slope  of  stress- 
strain  curve)  to  the  curves. 

5.  Procedure 


a.  Obtain  pull-test  data  for  particular  nylon  tape  from  the 
NAVAIRTESTFAC  Engineering  Department  Recovery  Division  (load  lb  - 
stretch  in.). 


b.  Convert  data  to  a  stress-percent  tape  elongation  curve. 

c.  Curve  fit  best  3rd-degree  curve  through  daca  to  obtain  equa¬ 
tion  that  represents  stress  as  a  function  of  strain. 
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d.  Differentiate  equation  to  obtain  a  new  equation  that  repre¬ 
sents  the  slope  of  the  stress-strain  curve  and  is  the  modulus  of  the  tape 
as  a  function  of  tape  elongation. 


E-4 


NATF-EN-1120 


APPENDIX  F  -  DETERMINATION  OF  TAPE- REEL  ACCELERATION 

1.  Tape-reel  acceleration  is  one  of  the  key  "motion  producing" 
parameters  of  the  program  (see  Simulated  Arrestment  Description,  Section 
V  of  main  text) .  The  importance  of  this  parameter  necessitates  an  ex¬ 
planation  of  the  equations  that  represent  it. 

2.  Tape-Reel  Acceleration:  The  equation  for  angular  acceleration 
of  the  tape  reel  is  obtained  from  the  relationship: 

ST  =  lot, 

where  ST  =  sum  of  the  torques  acting  on  the  reel  and  rotor, 

I  =  polar  moment  of  inertia  of  the  rotating  reel, 

and  ot  -  angular  acceleration  of  the  reel  and  rotor. 

a.  The  ST  term  is  the  sum  of  the  torque  created  by  the  tape 
tension  force  applied  at  the  radius  of  the  outer  wrap  of  the  tape  (Ttt) 
and  the  retarding  torque  developed  in  the  water  brake  (Twb) 

or  ST  -  Tj;t  Twb 

(1)  Ttt  =  Tt  Rt 

where  Tt  =  purchase-tape  tension 

and  Rt  =  radius  of  outer  tape  wrap  on  the  reel. 

(2)  Twb  =  B-  wn 

where  B  =  hydrodynamic  brake  constant, 

w  =  angular  velocity  of  reel  and  rotor, 
and  n=  brake  exponent  (constant). 

The  two  constants,  B  and  n,  are  obtained  from  actual  test  data  of  a  par¬ 
ticular  type  gear.  When  the  value  of  angular  velocity  w  reaches  its 
maximum  value  for  a  test  event,  then  the  angular  acceleration  is  mini¬ 
mized.  Assuming  that  it  becomes  zero,  a  =  0  and  w  =  wm  (maximum  angular 
acceleration),  from  the  original  equation, 

ST  =  let  =  0 

so  that  ST  =  Tj.  Rt  +  B  wm  =  0, 

or  Tt  Rt  =  B  w  • 
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The  values  of  Tt  and  Rt  are  obtained  from  the  test  data  at  the  time  cor¬ 
responding  to  the  maximum  angular  velocity  Wm.  The  values  of  B  and,  n 
are  then  obtained  by  using  the  method  of  least  squares. 

b.  The  X  term  is  the  sum  of  the  tape-reel  (metal)  inertia  Im 
and  the  instantaneous  tape  inertia  It  I  *  Im  +  If 

(1)  The  value  of  tape-reel  inertia  1^  is  a  constant  value  and 
is  equal  to  114  slug  -  ft  for  the  E-28  arresting  gear. 

(2)  The  value  of  the  tape  inertia  It  varies  as  the  tape  is 
pulled  off  of  the  reel  and  is  obtained  by  the  expression 

It  -  ^  (rh2  +  Rt2) 

where  M  =  instantaneous  tape  weight  /  g, 

Rjj  =  reel  hub  radius, 

and  Rt  *  instantaneous  radius  of  outer  tape  wrap. 

c.  Upon  rearrangement  of  the  described  terms  in  the  original 
equation,  the  expression  for  tape-reel  acceleration,  is: 

_  £T  _  Ttt  +  Twb  _  Tt  Rt  +  Bw11 
*  ‘  1  ’  IfI  +  Ic  '  Im  +  |m  (Rh2  +  Rt2)  ' 


> 
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This  mathematical  model  is  designed  to  provide  predicted  dynamic  performance 
data  of  shorebased  rotary-hydraulic- type  aircraft  arresting  gears.  A  Navy  Model 
E-28  arresting  gear  is  used  for  specific  comparison  between  predicted  results  of 
the  computer  solution  and  actual  test  results.  The  simulation  of  an  arrestment 
of  a  vehicle  under  a  particular  set  of  conditions  is  accomplished  by  putting  in¬ 
formation  (data)  into  the  computer.  The  input  data  specifies  values  for  the  in¬ 
stallation  geometry  and  mechanical  properties  of  the  arresting  system  and  the 
test  vehicle.  Predicted  dynamic  values  of  forces  and  motions  of  the  test  vehicle, 
purchase  system,  and  tape  reel  are  printed  out  versus  time  at  a  predetermined 
incremental  time. 


This  report  is  a  phase  report  on  the  development  of  the  model  and  contains 
the  early  analytical  design  approaches,  the  most  current  analytical  approach 
with  the  computer  program,  and  instructions  for  execution  of  the  computer 
program.  A 
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